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abstract 


trSce^flt:  install ed^^the^three  si  reP°rted'  Varlab'-“P-itancc 

broad-band  output.  aTcm/ \£rZlrlg^ T^/sTr^Thu", *"  ‘ 

second  325  deg  azimuth  strain  seismometer  will  he  „c-a  *«  K p  ‘ 

seismometer  in  July.  Magnets  and  coils  on  *he  id  •  d  ?°  construct  a  vertical 

saoo til ingCandb differentiating 

demonstrate  t'haTJn pT'lTL ITT^V  "“"-P***  --rdlgs 

nr  a 7  .  .  ..  ’  i-8«  ^n,  and  Sg  travel  times  from  Rangely,  Colorado  to 

s:»dbetwr  sssrr*  a”d  p.ia’Lsr. 

strains  a  «.  • ,  .n0  se  1S  reviewed.  Theoretical  equations  are  developed  for 

50  sec  ri  3  half"sPace  from  a  P^ne  wave  pressure  fluctuation.  At 

-a--  orizontal  inertial  seismometer  operating  at  TOOK  fat  25  cPri 

“rfoT^oVirTrir  45  •°°°  -  f°r 

ULl  u  27’900*000-  A  hypothesis  is  presented  that  at  all  ha^i  r0ck 
seismograph  stations  non-linear  rigid  block  tilting  occurs  as  a  result  of 

iseSSm  -IWn?  3nd  ^y,l6igh  WaVes-  Evidence  that  supports  thi  s  “hypothesis 
is.  (1)  TTiere  is  a  slight  relative  amplitude  shift  with  frequency  and  Tlo 

himr  PhdT  bet“een  11,6  pa«llel  st™n  sei,™^2S  lb. 

■ISSTSlSSri-.  reSP°nd  t0  »re  ti  it  than  is  assorted  . 
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DEVELOPMENT  OF  LP  WAVE  DISCRIMINATION 
CAPABILITY  USING  LP  STRAIN  INSTRUMENTS 
Quarterly  Report  No.  8,  Project  VT/8706 


1.  INTRODUCTION 


dl5(:ussesi.the  P^gress  during  April,  May,  and  June,  1970  in  the 
installation  and  preliminary  operation  of  long-period  (LP)  strain  and  inertial 
seismographs  with  equivalent  magnification  and  response  characteristics  The 

"  ‘  6ntS  are  to  develop  techniques  for  discrimination  of  LP 

seismic  waves.  The  major,  effort  on  e*ch  task  of  the  Statement  of  Work  is 

v  ^SSG  ln  SeParate  sections.  Thi£  report  is  to  apprise  the  Project  Office 
of  the  current  status  of  Project  VELA  T/8706.  It  is  submitted  in  Compliance 

eg^-Om611"  °f  ^  C°ntraCt  °ata  List,  Contr^t  ^3657- 


2 •  DEVELOP  DESIGN  SPECIFICATIONS.  TASK  a(l) 


This  task  has  been  completed. 


3<  DETERMINE  THE  MOST  EFFECTIVE  TECHNIQUE,  TASK  a(21 


This  task  has  been  completed. 


4-  DESIGN,  FABRICATE,  AND  TEST  LABORATORY  MODELS.  TASK  b(l) 

r 

This  task  has  been  completed. 


5. 


DEVELOP  A  FINAL  ENGINEERING  MODEL  DESIGN,  TASK  b(2) 


’  ^HlC^  "T  PTe^lously  reported  complete,  >as  reopened  when  it  became 
*  tha*  a  dr*ft  Pr°blem  in  the  Optical  Displacement  Transducer,  Model 
^'°v  -r^Uir  extensive  ^sting  and  possible  redesign.  Since  a  recently 
designed  Variable  Capacitance  Displacement  Transducer  (VCT)  had  been  success-  ' 
fully  operated  by  another  group,  it  was  decided  that  design  efforts  could  be 
most  profitably  directed  toward  an  adaptation  of  the  VCT  for  use  at  Queen 
Creek,  Arizona,  Seismological  Station  (QC-AZ) . 
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5.1  OPTICAL  DISPLACED  TRANSDUCER,  MODEL  32770 

SH  SENSES 

ZZ^-ZlTnlT^Z  £fM  W’ 

QC-AZ  instrument  nomenclature,  recorder A  mill \  Jn  seismom^ter-  CTTie 

ssssi-s rr^  „ 

ttest™?,1«ureabt1he"iriftnTirt0n”K"5'-  °ptical  elen,e"ts 

X.  and  exceeded  the  ^amic 

=“-=Si£S?:i 

Although  the  drift  was  obviously  an  effPrt  nf  +  J,  ^  1  M1  er  fgain- 
were  unsuccessful  *  i  ninDointino  t-i,  emperature  change,  field  tests 

the  optical  tra^sducerr  aLd8oned.SOUrCe‘  ^  de-loP-nt  of 


5.2 


VARIABLE-CAPACITANCE  DISPLACEMENT  TRANSDUCER 


about  the  S  SethLl ^A^^^rned  TT  “ 

™*-t  ^aJCT  described  by  Shc^cS**)!  r 

March,  and  April  o/Jhis  yearC1hl!enefLlddr  Pr°ject, VT/8703  du™g  February, 

were' made:  U  “*  ^ 

1.  Only  one  variable  capacitance  detector  was  used  for  the  transducer; 
QC-AZ  power°subsys tem^  CirCU‘tS  ""  “  be  “"Fatible  with  the  existing 

3.  Remotely  operated  coarse  and  fine  adjusting  assemblies  were  omitted; 

A.  A  10:1  parallelogram-lever  type  motion  reducer  driven  hv  a  7  „m  „„ 
iv  mlc^t'eT*81  ““  SUbStitUted  f°r  *ha  >  reducer  driven  by  a  W  Z, 

^si™5UcannbedhnifJ'  alti;OU8h  laCMn8  S°me  °f  the  features  »f  the  original 
aes 1 gn ,  can  be  built  quicker  at  a  considerably  lower  cost. 


-2- 


TR  70-29 


Preliminary  specifi cations  of  the 
Sensitivity  (detector) 

Noise  level  fdetector) 


transducer  are  as  follows: 


40  mV  Der  micron  displacement 


50  mV  p-p  in  the  period  range  of 
1000  to  1  sec 


Linear  range  of  operation 
Dynamic  range  (detector) 
Long-term  drift 
Frequency  response 
Cut  off  rate 
Min.  detectable  strain 
Power:  Voltage 

Power  • 


±10  Mm  minimum  for  ±1  percent  linearity 
70  dB,  minimum 
150  mV/°C 

Flat  in  strain  from  dc;  6  dB  down  @7.5  Hz 
12  dB/octave 

3  x  10-H,  dc  to  5  Hz  (40  m  rod  length) 

+19  to  +2 5  Vdc,  -19  to  -23  Vdc,  center 
common 

2  watts,  maximum 


S3-.SS-.  ~  b‘£v™  =.:rr 

between  the  other  plates  caDacitL.P  nr  k  ^  W  th  the  center  plate  midway 
In  the  detector,  Z1  and  Z2  aL  in^r*?^  H  caPacltors  is  nominally  100  pF. 

connected  as  a  free-running  multivibrator  “ArbllScr'ttT1  ‘T*""5 
outputs  of  Z1  and  Z2  are  tumotr  i  '  r  balance,  the  complementary 

If  the  center  plate  is  display? “^Ch  contain  »»  *=  component, 
capacitor  goes  down  while  the  orter  caPaaata"ce  of  one  side  of  the  dual 
displacement.  Proportional  chan^  ^apacitance  80es  UP  m  proportion  to  the 
used  to  control  multivibrator  symmetry6  time  C°nStant  circuits 

nonsymmetrical  square  wave  at  aDDroxfLr  i  e  ’"ultlvlb^ator  output  is  then  a 

balance.  These  nonsymmetrical  oStp^wre  then  fMt  that  at 

components  are  applied  to  the  positive  anrl  hen  f*ltered  and  the  resulting  dc 
The  output  of  Z3  is  a  dc  vnVr*  negative  inputs  of  amplifier  Z3. 

u  put  or  Lt  is  a  dc  voltage  proportional  to  displacement. 

completedf^No  Jurih8er1es?^SSe:%Ch\"itl'k“m7“t:df°^  ^  VCT  «“ 


6-  DESIGN  A  FIELD  TEST  INSTALLATION,  TASK  b(3) 


This  task  has  been 


completed. 
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7.  DESIGN,  FABRICATE,  AND  INSTALL  THE  FIELD 
TEST  INSTALLATION,  TASK  c(l) _ 


7.1  FABRICATE  INSTRUMENTATION 

When  design  work  on  the  VCT  was  completed  in  early  June,  all  necessary  electronic 
components  were  ordered  and  fabrication  of  the  mechanical  parts  was  begun. 

The  three  transducers  are  scheduled  for  completion  in  early  July,  after  which 
no  further  work  under  the  sub-task  is  planned. 


7.2  INSTALLATION  OF  THE  STRAIN/ INERTIAL  COMPLEX 

7.2.1  Additional  Equipment  Installed 

During  this  reporting  period,  several  items  of  equipment  were  added  to  the 
QC-AZ  system.  The  equipment  was  furnished  to  Project  VT/8706  as  Government 
Furnished  Equipment  (GFE) . 

The  first  major  addition  was  a  complete  microbarograph  system  for  monitoring 
outside  air  pressure  changes.  The  system  was  installed  on  8  April  1970  and 
has  been  recording  on  film  and  tape  as  ML2.  The  microbarograph  transducer 
senses  pressure  changes  with  a  single  port  and  the  overall  response  of  ML2  is 
flat  to  rate-of-change  of  pressure  from  about  47  to  200  seconds.  Calibration 
of  the  two  microbarograph  systems  outside  and  inside  the  mine  shows  that  the 
mine  seal  provides  more  than  40  dB  attenuation  of  outside  air  pressure  changes. 

On  16  May  1970  routine  operation  of  a  second  Develocorder,  Modwl  4000,  was 
begun.  Prior  to  putting  it  in  operation,  the  unit  was  dismantled  for  thorough 
(.leaning,  general  maintenance,  and  updating  with  a  light-slash  time  line  marker. 
Recordings  of  various  on-line  sums  and  differences  of  the  long-period  strain 
and  pendulum  seismographs  provide  four  beams  steered  toward  azimuths  of  325, 

55,  145,  and  235  degrees.  Data  are  being  used  to  demonstrate  the  on-line 
capability  of  the  strain-inertial  complex  to  enhance  surface  waves  and  to  aid 
in  wave  discrimination  and  wave  type  identification.  Signal  isolation  and 
gain  control  for  the  various  sums  is  being  provided  by  operational  amplifiers 
in  an  Analog  Computer,  EAI,  Inc.,  Model  TR-10.  The  computer  and  the  data 
control  equipment  for  the  second  Develocorder  required  addition  of  a  fifth 
equipment  rack  in  the  primary  recording  van. 

7.2.2  Future  Plans 

Data  collected  from  the  two  side-by-side  325  deg  azimuth  strain  seismometers 
are  considered  adequate  for  the  purposes  of  the  test.  Therefore,  the  temporary 
S325  seismometer  will  be  removed  and  its  components  installed  for  vertical 
operation  in  early  July.  Also,  it  is  expected  that  the  broad-band  transducers 
can  be  completed  and  installed  in  late  July.  At  that  time,  the  installation 
phase  at  QC-AZ  will  be  complete. 
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8.  CONDUCT  NOISE  AND  STABILITY  TESTS,  TASK  c(2) 


Noise  and  stability  tests  are  being  conducted  as  part  of  the  system  operation. 
Results  of  the  noise  and  stability  tests  are  given  under  tasks  d(l)  and  d(2) . 


9.  CONDUCT  PRELIMINARY  OPERATION  AND  EVALUATE 
INSTRUMENT  PERFORMANCE,  TASK  c(3) _ 


Preliminary  operation  has  been  completed  and  evaluation  of  instruments  will  be 
a  continuing  part  of  the  system  operation.  Instrument  evaluation  is  reported 
under  tasks  d(l)  and  d(2)  . 


10.  OPERATE  THE  STRAIN- INERTIAL  SYSTEM,  TASK  d(l) 


Routine  operation  of  the  strain-inertial  complex  was  continued  on  a  7-day-a- 
week  basis.  At  the  end  of  this  reporting  period,  21  separate  seismograph 
channels  and  5  environmental  data  channels  were  being  routinely  recorded  on 
two  16-mm  film  recorders,  one  35-mm  film  recorder,  and  three  14-channel  FM 
tape  recorders.  Data  were  interrupted  several  times  for  equipment  installation, 
modification,  ard  maintenance. 

10.1  EQUIPMENT  MODIFICATIONS 

10.1.1  Long-Period  Inertial  Seismometers 

In  early  April,  the  amplitude  responses  of  the  long-period  (LP)  inertial 
seismographs  indicated  a  serious  variation  from  the  theoretical  curve.  Since 
the  strain  seismographs  were  built  to  match  the  theoretical  Advanced  Long- 
Period  Seismograph  (ALPS)  response,  the  different  inertial  response  resulted 
in  a  degradation  of  on-line  sums  and  differences  between  strain  and  inertial 
seismographs.  The  problem  was  traced  to  the  seismometers  which  were  equipped 
with  more  efficient  transducers  than  the  standard  design.  These  prototype 
transducers  were  installed  in  another  program  prior  to  the  transfer  of  the 
seismometers  to  Project  VT/8706.  The  higher  magnet  charge  and  greater  length 
of  wire  in  the  coils  resulted  in  a  seismometer  with  a  much  higher  critical 
damping  resistance  (CDR)  than  normal.  For  example,  the  two  horizontal  seismo¬ 
meters  had  a. CDR  of  about  2500  0  as  compared  to  the  standard  value  of  1510  ft. 
This  high  CDR  resulted  in  overdamped  seismometers  when  other  system  parameters 
were  set  to  nominal  values.  In  the  ALPS  system  where  the  seismometer  and 
galvanometer  are  direct  coupled,  seismometer  and  galvanometer  damping  plus  the 
coupling  coefficient  are  interrelated  so  that  one  parameter  cannot  be  changed 
without  affecting  the  other.  It  was  therefore  necessary  to  modify  the  trans¬ 
ducers  to  obtain  standard  values  of  coil  wire  length,  coil  resistance,  and 
magnet  flux  density. 
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New  coils  with  less  wire  were  wound  and  installed  on  the  seismometers  beginning 
on  21  May.  These  new  coils  incorporated  a  recent  improvement  in  LP  seismometer 
coils  developed  especially  for  high-gain  operation.  This  improvement  balances 
the  dimagnetic  and  paramagnetic  properties  of  the  coil/coil-form  assembly  so 
that  non-linearities  in  restoring  force  are  minimized.  In  addition  to  the 
installation  of  new  coils,  magnet  charge  was  adjusted  on  the  seismometers. 

The  vertical  seismometers  presented  more  of  a  problem  than  the  horizontals. 

Since  the  vertical  instrument  requires  a  higher  flux  density  than  the  hori¬ 
zontals  (0.195T  and  respectively)  and  the  prototype  magnets  could  not 

be  charged  to  a  0.195T  flux  density,  standard  magnets  were  used  as  replacements. 
The  modifications  to  the  vertical  seismometer  were  completed  on  18  June.  A 
final  (fine)  matching  between  the  strain  and  inertial  seismographs  will  be 
made  in  the  near  future. 

10.1.2  Strain  Seismographs 

When  using  broad-band,  solid-state  amplifiers  with  velocity  transducers,  high 
frequency  signals  outside  the  normal  passband  can  very  easily  saturate  the 
amplifiers.  In  the  strain  system,  this  problem  was  anticipated  and  the  flag 
output  was  incorporated  in  the  strain  filter  unit  to  sense  this  excess  voltage. 
After  a  spare  20-channel  galvanometer  bank  became  available  to  replace  the 
galvanometer  bank  in  Develocorder  No.  1,  recording  of  the  flag  trace  as  a  20th 
channel  was  begun  on  27  April.  It  was  immediately  obvious  that  many  of  the 
large  spikes  recorded  on  LP  strain  channels  were  the  result  of  Ithaco  amplifier 
clipping  associated  with  high  amplitude  short-period  signals.  To  minimize 
this  preamplifier  clipping  problem,  capacitors  were  put  across  the  transducer 
coil  output  leads  on  11  May.  The  capacitors  and  coil  resistances  form  single 
pole  low-pass  filters  with  a  cut-off  comer  at  0.8  Hz.  These  filters  reduce 
10  Hz  signals  by  about  20  dB,  but  do  not  affect  periods  long  compared  to  1.25 
sec  (0.8  Hz).  The  filters  have  been  very  effective  in  reducing  the  frequency 
of  occurrence  of  saturation.  It  will  be  necessary  to  modify  the  short-period 
section  of  the  Strain  Filter,  Model  33350  to  achieve  matched  responses  between 
short-period  strain  and  inertial  seismographs.  The  modification  consists  of 
changing  the  normal  0.8  Hz  narrow  band-pass  filter  to  a  0.8  Hz  high-pass  filter. 
This  filter,  together  with  the  low-pass  coil -capacitor  input  filter,  will  form 
a  composite,  critically  damped  band-pass  filter.  It  will  be  necessary  to 
adjust  the  existing  circuits  to  change  the  damping  of  the  short-period  inertial 
seismometers  from  0.69  to  1.0  of  critical  to  achieve  matched  short-period 
responses . 

In  another  modification,  the  responses  of  the  ultra- long-period  (ULP)  strain 
seismographs  were  changed  to  reduce  the  signal  level  at  the  recorders  for 
periods  longer  than  about  300  seconds.  The  original  UPL2  response,  which  was 
flat  to  ground  displacements  from  travelling  waves  at  periods  from  30  to  1000 
seconds,  allowed  recording  of  high  level,  very  long-period  noise  which  limited 
the  usefulness  of  the  ULP  recordings.  The  new  response,  designated  ULP3,  was 
achieved  by  substituting  283  sec  low-pass  filter  cards  for  the  1000  sec  filter 
cards,  and  is  flat  to  ground  displacements  from  travelling  waves  at  periods 
from  30  to  283  seconds.  The  283  sec  low-pass  filter  is  used  in  the  ULP  1 
response,  so  the  cards  for  this  change  were  on  hand.  In  addition  to  the  283 
sec  filter  change,  the  30  sec  low-pass  filter  cards  were  modified  to  provide 
an  additional  14  dB  of  circuit  gain.  This  was  necessary  to  obtain  sufficient 
background  voltage  for  magnetic  tape  recordings. 
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10.2  MAINTENANCE 

During  April,  considerable  time  was  devoted  to  locating  the  cause  of  excessive 
n^se  on  the  LP  inertial  seismographs,  especially  P55L.  The  noise  was  found 
to  be  a  combination  of  ground  loops  and  60  Hz  voltage  pickup  m  the  “"amplified 
data  circuits.  Ground  loop  problems  at  QC-AZ  are  mc^ased  by  higlr  ground 
resistivity  which  results  in  a  60  Hz  potential  difference  of  about  10  \  P  p 
between  the  seismometer  room  and  the  adit  of  the  mine.  This  problem  was 
minimized  by  careful  grounding  techniques.  The  60  Hz  voltage  pickup  was  foun 
to  "cur  i7 the  twisted-pair.  No.  12  AWG,  solid  copper  data  Unes  whrch  extend 
some  460  feet  from  the  seismometers  to  the  phototube  amplifiers  (.PlA  sj  . 

Shielded  data  lines  solved  the  problem. 

I  iehtning  damaged  the  Ithaco  amplifiers  in  the  strain  system  on  27  April  and 
on  24  June.  In  both  cases,  shorted  transistors  in  the  amplifier  output  stages 
indicated  that  overvoltages  were  occurring  on  the  field  lines  between  th 
amplifier  and  the  vans.  The  lightning  protection  circuits  consist  of  gas 
discharge  diodes  and  high  surge  semiconductor  diodes  across  the  lines,  and  laJ~ge 
inductors  in  series  with  the  lines.  These  circuits  are  identical  with  circuits 
being  successfully  used  on  other  programs  but  no  reason  has  been  foun 
their  failure  to  adequately  protect  the  QC-AZ  circuits. 

During  this  period,  the  mine  seal  was  improved  by  eliminating  leaks  around  the 
second  and  third  ship’s  doors  in  the  mine.  Larger  leaks,  especially  aroun 

the  masonry  holding  the  doors,  were  sealed  with  caulkJn«  ^  '-ere 

leaks  were  found  in  the  many  cracks  in  the  rock  around  the  doors  These  were 

sealed  by  painting  a  relatively  large  area  around  both  sides  of  the  doors 
a  Portland  cement  slurry.  The  time  constant  of  the  instrument  room  behin 
the  third  door  is  now  about  4  hours  and  the  time  constant  of  the  area  between 

Te  second  and  third  doors  is  about  40  minutes.  This  is  considered  more  than 

adequate  for  successful  operation. 

10.3  AMPLITUDE  AND  PHASE  RESPONSE 

r;;:s,= 

an7p!o«ed  atSa7ateran«,e.  Phase  responses 
are  obtained  by  adjusting  phase  of  a  calibration 

zero'ph  as  (^between  ST  Z&ZT*.  sei  saograph 

ourPu?  After  nu^eroL  cycles  a?e  observed,  phase  adjusted,  and  reobserved, 
the  phase  difference  is  read  from  the  dial  and  written  own. 

and^ralelrt  square  ^/between  a  Seoreti cruise  from  a  seismograph  with 
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it  required  less  lost  data  than  the  first  method  and  because  it  produced  a  more 
satisfactory  solution  than  the  third  method. 

To  obtain  the  seismograph  transfer  function  from  a  random  noise  input,  a 
sufficiently  long  data  sample  must  be  obtained  to  provide  the  desired  statisti¬ 
cal  stability  of  the  cross  power  spectral  density  between  the  input  and  the 
output.  A  data  sample  at  least  32  min  duration  would  be  necessary  for  periods 
up  to  100  sec  and  at  least  6.7  hr  duration  would  be  necessary  for  periods  up 
to  2000  sec.  An  attempt  was  made  to  use  the  white  noise  from  a  resistor  on 
the  input  of  a  broad-band  dc  amplifier  driving  the  calibration  coil  of  an 
inertial  LP  seismometer.  With  the  broad-band  system,  60  Hz  noise  dominated  the 
calibration  coil  input.  When  it  became  evident  that  extreme  grounding  pre¬ 
cautions  would  be  necessary,  the  experiment  was  deferred  to  the  other  two 
methods  and  to  other  more  pressing  installation  problems. 

The  three  inertial  P(Z,  325,  55)L  seismographs  were  attenuated  to  improve  tne 
signal-to-noise  ratio  and  several  steps  of  current  were  applied  to  each  of 
the  three  calibration  coils  separately.  Figure  2  is  l  recording  made  from  the 
magnetic  tape.  It  shows  the  reduction  in  microseismic  background  by  the 
attenuation  and  the  first  current-on  and  current-off  and  second  current-on 
pulses  on  the  P55L  seismograph.  These  pulses  were  digitized  and  were  used  for 
evaluation  of  both  methods  two  and  three.  Both  methods  use  the  fast  Fourier 
transform  algorithm  of  Cooley  and  Tukey  to  perfornrAhe  required  calculations 
rapidly. 

The  second  method  is  straightforward,  but  contains  the  problems  associated  with 
numerical  differentiation.  The  seismogram  pulses  are  equivalent  to  the  output 
of  the  seismograph  for  a  step  of  acceleration  input.  (A  step  of  current  into 
the  calibration  coil  produces  a  step  of  force  F  on  the  inertial  mass  M  which 
is  accelerated  at  the  rate  of  F/M.)  The  desired  amplitude  and  phase  response 
is  for  an  input  of  an  impulse  of  displacement.  To  obtain  the  desired  transfer 
function  three  derivatives  are  taken.  Because  of  the  excellent  signal-to-noise 
ratio,  the  fi  -st  attempt  was  to  take  the  Fourier  transform  into  the  frequency 
"domain  and  take  a  triple  derivative  by  multiplying  by  (iu>)  .  The  first 
derivative  (figure  3a)  appeared  smooth  at  all  frequencies,  and  the  desired 
third  derivative  (figure  3b)  was  relatively  smooth  at  periods  from  512  down  to 
20  sec.  Below  20  sec  the  noise  caused  erratic  results.  A  seven-point  parabolic, 
least  squares,  smoothing  (Lanczos,  1956)  was  then  accomplished  before  each 
derivative  was  taken.  The  smoothed  amplitude  response  is  plotted  down  to  10 
sec  in  figure  4  along  with  the  experimentally  measured  amplitude  curve.  At 
18  sec  and  12  sec  microseisms  are  seen  to  contribute  noise  to  the  result;  how¬ 
ever,  when  matched  with  the  experimental  curve,  the  results  are  very  good  when 
the  microseismic  energy  is  allowed  for.  The  phase  response  curve  is  plotted 
in  figure  5.  A  five-point,  parabolic,  least-square,  smoothing  was  also  tried, 
but  the  true  data  are  sufficiently  smooth  that  the  seven-point  smoothing 
produced  better  results  for  very  little  additional  computation  time. 

The  third  method  used  to  obtain  the  seismograph  transfer  function  was  considered 
unsuccessful.  The  method  was  developed  by  Mitchell  and  Landisman  (1969)  for 
the  World  Wide  Standard  Seismograph  Network  (WWSSN)  recordings.  A  set  of 
seismograph  paramete’s  a tie  input  to  the  program.  A  theoretical  time  domain 
response  to  an  input  step  of  acceleration  is  calculated  using  numerical  inte¬ 
gration  which  is  inherently  more  stable  i  han  numerical  differentiation. 
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RELATIVE  DISPLACEMENT  AMPLITUDE  RELATIVE  ACCELERATION  AMPLITUDE 


10 


a.  P55L  response  to 

impulse  of  accelera¬ 
tion  (one  deriva¬ 
tive,  no  smoothing) 


PERIOD  (seconds) 


10 


b.  P55L  response  to 

impulse  of  displace¬ 
ment  (three  deriva¬ 
tives,  no  smoothing) 


Figure  3.  Amplitude  response  of  P55L  calculated  from  the  pulse  in  figure  2 
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KEIATIVE  DISPLACMENT  AMPLITUDE 


Figure  4. 
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Figure  5.  Phase  response  of  P55L  calculated  from  the  pulse  in  figure 
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A  set  of  partial  derivatives  of  the  time  domain  pulse  to  changes  in  seismometer 
and  galvanometer  free  period  and  damping  is  calculated.  A  least  squares  adjust 
ment  increment  in  these  four  parameters  is  calculated.  The  parameters  are 
adjusted  and  the  program  iterates  until  the  residual  sum  of  squares  between  the 
theoretical  and  the  actual  time  series  is  below  a  predetermined  acceptable 
level.  The  WWSSN  LP  seismograph  has  a  seismometer,  an  attenuation  circuit  a 
galvanometer,  and  a  photographic  recorder  with  a  light  beam  reflected  off  the 
mirror  in  the  galvanometer.  The  QC-A2  seismographs  contain  electronic  filters 
that  do  not  have  a  counterpart  in  the  WWSSN.  It  was  necessary  to  add  the 
transfer  function  of  these  filters  to  the  program  so  that  the  proper  theoreti¬ 
cal  response  could  be  calculated.  The  filters  are  PTA  plug-in  filters,  Model 
No.  6824-14,  without  a  6  sec  notch,  and  Model  No.  6824-15,  with  the  notch. 

The  parameters  of  the  actual  filters  vary  slightly  from  one  to  another,  and 
the  use  of  slightly  incorrect  comer  frequencies  may  have  affected  th?  results 
some,  but  more  serious  problems  Were  encountered.  Since  Mitchell  and  Landisman 
had  such  good  success,  the  program  was  first  used  with  real  data.  After 
completely  unacceptable  results  were  obtained,  the  measured  system  parameters 
were  input  to  the  program  and  the  theoretical  time  domain  pulse  was  output  to 
use  as  an  artificial  signal  for  a  test  case.  With  this  artificial  signal, 
and  parameters  within  5  percent  of  the  true  parameters  input  to  the  program 
it  did  not  arrive  at  a  correct  solution  after  25  iterations,  although  it  was 
slowly  incrementing  in  the  right  direction.  Because  the  program  did  not 
converge  rapidly,  because  it  was  very  sensitive  to  the  pulse  start  time  and 
thus  signal-to-noise  ratio,  and  because  of  the  longer  computation  time,  this 
method  was  abandoned  in  favor  of  the  second  more  straight-forward  method. 


10.4  FUTURE  PLANS 

Out  of  necessity,  operations  have  been  continued  throughout  parts  of  the 
installation  and  testing  phases.  As  a  result,  data  quality  has  suffered 
somewhat,  although  collected  data  are  adequate  for  the  most  part.  In  the  next 
quarter,  work  on  installation  and  testing  tasks  is  expected  to  decrease. 

More  emphasis  can  then  be  placed  on  the  collection  of  consistent,  well  cali¬ 
brated,  high  quality  data. 


11.  DEVELOP  METHODS  OF  WAVE  DISCRIMINATION,  TASK  d(2) 


Methods  of  wave  discrimination  and  methods  of  verifying  proper  instrument  ' 
operation  are  being  developed  and  pursued.  In  this  report  the  preliminary 
results  obtained  with  some  of  these  methods  are  presented.  An  example  of 
the  enhancement  of  LP  seismic  waves  from  two  interfering  events  is  demon¬ 
strated.  An  illustration  is  presented  of  the  enhancement  of  the  arrivals  from 
an  explosion  and  of  the  arrivals  from  an  earthquake  in  the  coda  of  the  Lg  wave 
from  the  explosion.  Examples  of  detection  of  surface  waves  from  small  events 
including  an  =  4.2  at  an  epicentral  distance  of  100  deg  are  illustrated. 
Results  of  an  experiment  to  compare  parallel  horizontal  strain  seismometers 
are  discussed.  And,  linear  and  non-linear  tilt  noises  on  the  horizontal 
inertial  seismographs  are  analyzed. 
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11.1 


ENHANCEMENT  of  signals  from  interfering  events 


arrival  ti^ J^f^any^the^l^tirw41  °PP°Site  a2imuths  from  QC-AZ.  The 

<ujd  Geodetic  Survey  (USCSGS)  Preliminawltete™!13^'  United  States  Coast 
on  the  two  events  are:  "unary  Determination  of  Epicenter  (PDE)  data 

Bvent  1:  0  -  1003:42.2,  S5.0S,  2B.3W,  South  Sandwich  lslands  Region 

J  '  7°'  %  =  6'°-  4  =  113  5  *8'  a2i"U'h  *“*«  »  epicenter  .  J42,  de, 

0  -  2030:54.7,  SI. SN,  178. SW,  Andreanof  Islands,  Aleutian  Is 
h  =  48  - 


% 


5'7,  i  ‘  50-6  deS-  a2i”Uth  nation  to  epicenter  .  312.0  deg. 


p  •  -  cpn-enrer  =  312.0  de 

circle  paths'^  PD8*?"“,™>t‘™  ttreugho^t ' tt’e^ep'ort f"?.  the  PDE  d«a,  but  are 
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Figure  6.  Great  circle  paths  from  two  earthquakes  to  QC-AZ 
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Table  1.  Relative  response  of  strain/inertial 
_ complex  to  two  earthquakes _ 


LOCATION 

DISTANCE 

AZIMUTH  (STATI ON¬ 
TO  -EPICENTER) 


RESPONSE  IN  325  DEG 
AZIMUTH 
RAYLEIGH  WAVES 
LOVE  WAVES 

RESPONSE  IN  145  DEG 
AZIMUTH 
RAYLEIGH  WAVES 
LOVE  WAVES 

RESPONSE  IN  55  DEG 
AZIMUTH 
RAYELIGH  WAVES 
LOVE  WAVES 

RESPONSE  IN  235  DEG 
AZIMUTH 
RAYLEIGH  WAVES 
LOVE  WAVES 


EVENT  1 

SANDWICH  ISLANDS 

113.5  deg 

142.5  deg 

6.0 


0 

0 


2.0 

0 


0 

0.95 


0 

1.05 


EVENT  2 

ALEUTIAN  ISLANDS 
50.6  deg 
312.0  deg 

5.4 

1.9 

0 

0 

0 


0 

1.17 


0.25 

0.75 
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Figure  8.  Magnetic  tape  playback  of  arrivals  from  ■  t  * 

1:  0  =  2003:42.2,  south  Sandwich ^slL're!  "J"8  eVentS*  Event 

6.0,  A  -  113.5  deg,  azimuth  -  14?  c  a  S10n»  h  =  70  km,  mb  = 
Andreanof f  Island*  h  =  48  ^  =4^7  d?g;  cVent  2:  0  =  2°30:54.7 

de8  •  ,  A  =  50.6  deg,  azimuth  =  312.0 
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The  arrival  times  of  many  earthquake  phases  are  enhanced  considerably  with  wave 
identification  clarified  and  in  many  cases  the  time  breaks  are  sharpened  for 
accurate  timing.  Earthquake  phases  are  identified  in  the  figure  with  the 
phases  from  event  1  on  the  upper  line  and  from  event  2  on  the  lower  line. 

All  of  the  arrivals  from  event  2  are  smaller  than  those  from  event  1.  Several 
arrivals  are  particularly  interesting.  The  time  break  is  sharpened  for  PS 
and  PPS  on  S-P325L  at  2032:37  and  2033:23.  An  unidentified  phase  arrival  is 
detected  on  the  PZL  and  S-P325L  at  2034:08.  The  SSS  and  SKKKS  arrivals  are 
detected  on  S-P325L  and  can  be  timed  at  2043:24  and  2044:18.  The  arrival  of 
the  LQ  Love  wave  from  event  1  is  enhanced  on  the  S-P55L  at  about  2053.  The 
possible  arrival  of  LR1  Rayleigh  wave  from  event  2  at  about  2053  is  indicated 
by  possible  signals  on  S+P325L,  S+P55L  (?) ,  PZL,  and  S325L+S55L.  The  LR1 
Rayleigh  wave  from  event  1  at  2058  is  enhanced  on  the  S-P325L  and  cancelled 
on  the  S+P325L  even  with  the  magnetic  tape  recorder  clipping.  Many  body  phases 
are  detected  as  arriving  after  the  times  for  presently  tabled  phases.  These 
phases  are  probably  either  (a)  multiple  core  reflections,  (b)  second  reflected 
phases  arriving  by  the  dong  arc,  and/or  (c)  thrice  reflected  phases  arriving 
by  the  short  or  long  arc.  ^ihp  of  these  body  phase  arrivals  are: 


Suspected  type 
of  arrival 

Nearest 

minute 

Channel  detecting 
phase 

P 

2101 

PZL,  S325L,  S-P325L? 

S 

2101 

S325L-S55L,  S-P325L? 

P  or  S? 

2103 

PZL,  S325L+S55L 
S325L-S55L,  S-P325L 

P 

2104 

S-P325L,  S-P55L 
S325L+S55L ,  PZL 

P  or  S? 

2105 

PZL,  S-P325L ,  S+P55L 
S325L+S55L,  S325L-S55L 

P  or  S? 

2108 

S-P55L,  S32sL+S55L,  PZL 

P  or  S? 

2109 

PZL,  S-P325L ,  S-P55L 
S325L-S55L 

11.2  ENCHANCEMENT  OF  SHORT-PERIOD  WAVES 

Although  the  primary  objective  of  Project  VT/8706  is  to  develop  methods  of 
long-period  wave  discrimination,  magnetic  tape  recordings  are  being  made  of 
short-period  (SP)  inertial  and  strain  seismographs.  The  SP  strain  and  inertial 
seismographs  can  be  added  and  subtracted  for  the  same  azimuthal  response  as 
snown  in  figure  7.  Two  examples  of  enhancement  accomplished  with  the  short- 
period  systems  are  shown  in  figures  9  and  10.  At  the  time  of  these  recordings 
the  P55S  inertial  seismograph  was  inoperational .  Magnifications  for 
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Figure  12  illustrates  the  LQ1  Love  wave  from  a  magnitude  4.1  event  in  the 
Greenland  Sea  at  an  epicentral  distance  of  61.9  deg.  There  are  large  100  sec 
pressure  variations  causing  tilting  which  are  recorded  on  the  horizontal  inertial 
seismographs  and  ground  strain  recorded  on  the  strain  seismographs. 

Two  weeks  later  the  second  Develocorder  was  operational  and  on-line  sums  and 
differences  were  being  recorded.  Figure  13  is  the  recording  of  the  Rayleigh 
wave  train  from  a  magnitude  5.0  earthquake  near  New  Britain  Island  at  an 
epicentral  distance  of  98  deg  and  an  azimuth  of  270  deg.  This  magnitude  event 
from  the  south  Pacific  is  always  well  recorded.  Figure  14  is  a  reproduction 
from  camera  2  for  the  same  event.  At  the  time  of  this  recording,  the 
magnifications  had  not  been  determined  on  this  camera  and  the  strain  and 
pendulum  traces  were  not  equally  contributing  to  the  sums  and  differences. 
However,  the  Rayleigh  wave  is  seen  to  be  greatly  enhanced  on  the  S55L+P55L 
trace.  Both  the  S325L1  and  S325L2  were  not  operating  properly  during  this 
recording. 

Figure  15  illustrates  recording  of  a  Rayleigh  wave  train  from  a  magnitude  4.4 
earthquake  off  the  coast  of  Ecuador  at  an  epicentral  distance  of  42  deg  and 
an  azimuth  of  135  deg.  The  six  inertial  seismograph  traces  are  all  overlining 
and  the  two  325  deg  azimuth  strain  seismographs  are  responding  to  the  Rayleigh 
wave.  The  camera  2  recording  of  the  same  event  is  shown  in  figure  16.  The 
S325L2-P325L  difference  trace  is  showing  good  enhancement.  The  S55L  and  P55L 
seismographs  were  not  contributing  equally  to  the  sum  and  difference  and  the 
enhancement  does  not  show  as  well  on  these  traces. 

Figures  17  and  18  are  reproductions  of  recordings  from  cameras  1  and  2, 
respectively,  of  a  magnitude  4.2  earthquake  in  Turkey  at  an  epicentral  distance 
of  99.5  deg  and  an  azimuth  of  29.6  deg.  In  figure  17,  the  event  wohld  possibly 
not  be  detected  if  only  the  ALPS  response  inertial  seismographs  were  present. 
However,  the  30  sec  signals  are  seen  more  clearly  on  the  extended- long-period 
(XLP)  response  because  the  large  microseisms  present  at  this  time  have  been 
eliminated  by  the  XLP  response.  Even  so,  the  detection  of  this  event  would 
still  be  questionable  except  for  the  clear  enhancement  in  figure  18  on  the 
S325L1-P325L  and  on  the  S55L-P55L  difference  traces  and  the  clear  cancellation 
on  the  two  sum  traces.  Note  also  that  the  8  sec  microseisms  are  enhanced  on 
the  S55L+P55L  sum  and  generally  cancelled  on  the  difference  trace.  (The 
camera  drive  was  erratic  between  0610  and  0611.) 

Figure  19  is  a  copy  of  a  magnetic  tape  playout  from  a  recording  of  the  Rayleigh 
wave  from  a  magnitude  5.3  earthquake  in  the  Tonga  Island  Region  at  an  epi¬ 
central  distance  of  78.9  deg  and  an  azimuth  of  238  deg.  This  figure  illustrates 
the  90  deg  phase  difference,  in  accord  with  theory,  between  the  vertical 
inertial  seismograph  PZL  and  the  sum  of  the  two  orthogonal  strain  seismographs 
S325L1+S55L.  The  success  of  this  comparison  suggests  that  the  vertical  strain 
seismograph  will  be  a  valuable  addition  to  the  QC-AZ  instrumentation  complex. 
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Figure  10.  Enhancement  of  short -period 
arrivals  from  an  earthquake 
during  the  coda  of  an  earlier 
event  .-■..  1 
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s»."snss  N°- 2  f»  21  ^ 

0  ’  event ^•’CAN"27°  06’50'°"  N'  I16“  °<’S1.0"W,  Nevada  Test  Site  (NTS) 

h  ~  31 sm5  deg  8’  *  =  S’5  deg  “  613  kl"’  azilmith  station-to-epi center  = 

•me  vertical  strain  seismograph  is  formed  by  the  sum  of  the  perpendicular 
horttontal  strain  seismographs.  Veiy  clear  arrival  times  caTbTpicked  for 
Thpt;  nJtial  Pn  a,rlval  on  the  S-P325S  and  the  (S325S+S55S)  +PZS  traces 
fr  ?tra.CeS  sho"  ™ch  sharper  first  breaks  than  the  P325S  and  PZS  inertial 

arrival  of  ThTf.i  ^  Cl6arly  detected  between  the  Pn  and  Pg.  l^e 

r  .  e  ®  phases  at  1502:50  saturates  the  magnetic  tape  recorder 
erectromcs.  Several  distinct  P  and  several  distinct  S  arrivals  can  b°  seen 
in  the  coda  after  the  electronics  stopped  clipping.  Times  Were  picked  from 
a  high  speed  playout  and  the  travel  time  was  divided  by  the  distLce  to  obtain 

f0r  VariOUS  lU  b-»“"  ™  - 


Travel  time/distance 


Sg-Lgj 

Sg-Lg-, 


7.534  km/sec 
7.474  km/sec 
6.662  km/sec 
6.207  km/sec 
3.996  km/sec 
3.597  km/sec 
3.557  km/sec 


%  i 

During  the  Lg  coda  from  the  explosion,  arrivals  were  recorded  from  an  earth- 

neareR^elv  dlT Geolo«ical  SurveX  CUSGS)  close-in  network 
icnc.If?  ^noC  do  (See  flgure  10).  The  USC&GS  PDE  card  data  are:  0  =  . 

A  -  7  2 * dep4°  R043k'N’  1°8<>  u3*8™’  RangelX*  Colorado,  h  =  4  km,  =  4.6, 

A  -  7.2  deg  =  804  km,  azimuth  =  16.9  deg.  The  Pn  (or  P,)  phase  is  recorded 

with  an  impulsive  break  on  the  S-P325S,  S325S+S55S,  and^he  (S325S+S55S1  PZS 
traces  about  3  sec  before  the  Pg  phase.  Several  S  and  surface  wave  phased 
including  Sn  are  clearly  enhanced.  The  travel  time  divided  by  the  distance 
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Phase 


Travel  time/distance 


Pn(or  Pj) 

* 

6.101  km/sec 

Pg 

5.983  km/sec 

Sn  i 

5.655  km/sec 

Sn^ 

3.640  km/sec 

Sn3 

3.631  km/sec 

Sg-Lg! 

3.486  km/sec 

Sg-Lg2 

3.478  km/sec 

The  travel  times  and  velocities  between  Rangely,  Colorado,  and  QC-AZ  are 
considerably  slower  than  the  travel  times  from  NTS  to  QC-AZ.  To  the  authors' 
knowledge,  the  above  observations  are  the  first  reported  Sn  travel  times  in 
the  Basin  and  Range  Province. 


11.5  ENHANCEMENT  OF  SURFACE  WAVES  AND  DETECTION  OF  LOW  LEVEL  SIGNALS 

Several  examples  of  recordings  of  surface  waves  and  the  detection  of  low  level 
signals  are  illustrated  in  this  section..  These  examples  are  shown  as  a 
collection  of  individual  record  sections  and  do  not  represent  a  unified  study. 
During  much  of  the  time  covered  by  this  report,  installation  of  instruments 
was  the  primary  goal  and  perfectly  calibrated  recordings  with  equalized  gains 
and  correct  polarities  was  of  secondary  importance.  Accurate  location  of 
epicenters  is  not  available  for  several  weeks  after  the  events;  therefore, 
several  earlier  recordings  have  been  used  to  give  a  preliminary  demonstration 
of  the  QC-AZ  detection  threshold.  In  the  group  of  recordings  presented,  an 
improvement  in  record  quality  and  a  decrease  in  instrument  noise  can  be  seen 
as  time  progresses. 

Figure  11  is  a  recording  typical  of  the  many  earthquakes  that  occur  near 
Chiapas,  Mexico.  A  very  clear,  large  amplitude  slightly  dispersed  Love  wave 
on  the  P55L  .  id  P55X  precedes  the  larger  amplitude  Rayleigh  wave  on  the  PZL, 
,P325L  P.ZX  and  P325X.  The  S325L  had  been  set  at  too  high  a  magnification 
during  this  recording  period.  The  S55L  did  not  record  the  Rayleigh  wave 
because  it  is  almost  transverse  to  the  wave  motion  and  on  a  null  for  both  Love 
and  Rayleigh  waves.  The  large'  recorder  amplitude  at  which  this  magnitude  4.0 
event  from  an  epicentral  distance  of  25.1  deg  is  a  good  example  of  the 
capabilities  of  the  ALPS  response  inertial  system  which  has  a  magnification  of 
113  to  122K  in  the  X10  reproduction  in  figure  11.  The  ML2  microbarograph 
outside  the  mine  is  being  calibrated  with  a  65  ybar  input  at  112  sec  and  with 
a  9  ubar  input  at  4.9  sec  periods. 
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11.4  COMPARISONS  BETWEEN  THF  pad»,,c,  I  '/ 

WEEN  THE  pARALLEL  STRAIN.  SEISMOGRAPHS 

Til  1  Q  C  a  r-»4-  4  »  I 


#  - r^umvrno 

inis  section  describes  thp  i  *  I  f 

:SS"  t -S“-'sSSF^ 

-.p.visons  between  •*«»*  "‘t'SSSSf*, 

=rdrE“7^S,"« 

i?ah-i  r  SerieS  and  a  sma11  number  Jf  1^/°  not  coherence  squared  ) 

stability.  Time  samples  of  Mil*  5  la,gs  Were  USed  for  good  statist^!, 

bT  T'  “*«  f-XectL  beiween  if  Ia*S)  °r  32  0  U  wi  S 

6.7  hr  (with,  128  laps)  unm  .  d  256  sec  periods  and  timo  «r - 

-d  .  snaner.cSf^eea1i»Ct.ICUU,ed“ith  128  *' 

Sri'S-fiS  ^nth'dTsp^acemenfat^fseffor  to 

spectte,  s^feid  2s."  ■ i^r 

SS:^.^ £«><■£  ^ ^ V‘s^™ W  tbe 
calibration  period  of  25  or  ion  q  eslTec*  period  to  the  response  of  the  c  f*-10  ° ^ 
to  displacement  the  st-rai  i.  .  eC‘  ln  t^le  figures  in  which  str  •  -  y^tem  at  the 

r?- ^  issjUt^Z^  MffSWwS"' 
zzsz?  based z  brd  «»”«■> 

passband.  The  passband  tsTd  ^STa^TT*1  n™l^ed  toT?  Hz^' 

7.81  x  10-3  Hz,  1/256  =  T91  x *CtUal  emulations  was  either  1/128 
reader  wishes  to  normalize  to  a  Zl  °r  1/2560  =  3.91  x  l0-4  J 
with  several  published  m'  •  Pa?s6and  of  1  mHz  =  1/1000  H-7  fnT- 
divided  by  1000.  *  ««™oas»uc  spectra,  the  plotted  JaLs f  L^"”0" 
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Table  2. 

instrumentation 

systems  compared 

Signal 

condition 

System  1 

ti 

Sys  tenu  2 

Period  range  (sec) 

Figure  no. 

Low  wind 

J  S32SL1 

S325L2 

1-256 

21 

- 

S32SL1 

S325L2 

1-128 

22 

S325L1 

P325L 

1-256 

25 

S325L1 

Q 

P325L 

1-128 

°«9 

24 

r 

S325L1 

S55L,. 

1-256 

25 

High  wind 

S325L1 

S325L2 

J  1-256 

27 

S325L1 

S325L2 

20-2560 

28 

S325U1 

S325U2 

20-2560 

29  ; 

S325L1 

P325L 

1-256 

30 

S325L1 

S55L 

1-256 

31 

S325U1 

S55^ 

20-2560 

32 

Earthquake 

1 

1 

LR1  coda 

S325LL1 

f  S325LL2 

1-256 

33 

LQ3,  LR2,  and 
LR3 

S325LL1 

S325LL2 

1-256 

35 

LQ3,  LR2,  and 

1  LR3 

S325LL1 

S55LL 

1-256 

36 

LR1 

S325L.1 

S325L2 

1-256 

38 

LR1 

S325L1 

S325L2  * 

1-128 

39 

LR1 

S325L1 

P325L 

1-128 

40 

High  .Wind „ 

S325U1 

ML2 

20-2560 

42 

0  '■ 

S5SU 

» 

*• 

‘  °  „  *. 

ML2  ; 

20-2560 

•a 

43 

‘0 

,  t  ^  ■ 

'  *  *-}'  <»„  0 
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o*  .  ^ 
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Many  of  the  spectra  show  1  or  2  neaks  at  i  7  i  a  _ 

residual  powei/ noise  in  the  playback  and  3-‘2  S6C<  111656  peaks  are  the 

playback  £as  at  100  times  the  fiefd/  digitization  electronics.  Actual 

playback  relates  to  an  equivalent t.6  tiTtl  trUe  60  Hz  noise  at 

of  the  field  recorder.  °r  1,7  S6C  noise  in  the  time  scale 

*"e course,oi  the  “>«->«*•- »».  is 
Fourier  transform fFFTl  almrfih™  „r  rP  f  T  calculateJ  using  the  fast 

Geote*  CDC  ilOO^cSefra^e'S  M^IjVJYV"  the 

number  in  the  calculations  tUc  ^  t0  ^  below\the  largest 

FFT  algorithm  takes  differences  froiTalmos ^  ^  the  faCt  that  the 

used  employs  three  FFf  and  erratic  results  £umbef*;  1116  spectral  progran 

largest  spectral  number.  T  ^  ^ 

noise.  When  two  time  series  are  analv7»d  * Ultf  are  interpreted  as  computational 
program  using  the  FFT  thev  mult  ““^“d  together  with  the  cross  spectral 

>»5  to  106  nr|„geeofF„To  co^uSLIl^ote"  *°,2*k«  J  ”  "lth‘"  tha 

comparing  the  inertial  seismograph  and  th*  'ct  ®  first  cross  spectra  calculated 
since  the  strain  power  values  are  39  5  xk-l^smaH ismograph  showed  this  noise 
values.  In  subsequent  calculations  th/sJk  f  than  the  disPlacement 
this  numerical  problem.  Strain  ValU6S  Were  scaled  to  alleviate 

11-4.2  Low  Wind  Conditions 

SGleCted  ^r  the  comparison  when  the  two  para,,.!  st^in 
un°?*l  mann6r'  when  the  wind*  was  varying***'' 
motion  in  thTrfcorkgs  Thil  tilt  n°  appar6nt  eartMuake 

background  in  Mav.  Figure  20  is  a  nnrt'P  6  5ypi?al  of  the  QC-AZ  microseismic 
for  the  comparison.  Jhile  e  t nK”  °f  ^  tlme  int6rVal  that  W?s  used 
the  horizontal  inertial  seismL™ k u S615mo«rapl "ere  operating  normally, 

and  were  still  noisy,  not  having  settled  bOOk  t^n!,™  jagnets  and  coils  Ranged 
P325  seismometer  was  the  last  tn  htltt  °ack  to  normal  operation  yet.  Hie 
20,  the  S325L1  and  q*?ci  -j  orked  on  and  is  the  noisiest.  In  figure 

periods  2? similarities  at  the  shorter” 
motion  than  the  &325L1.  ’  S325L2  contains  more  longer  period 

s^L:„°*r^TurirSeSs%Je^^LVd  “  1  IT ^ 

end  1.764  x  10-24  strain2  on  the  Inn  f  10  stra1"2  »"  the  S325L1 

rsbands  ™  iTo'-nii^ 

a  lPt  sli;e  I:to^  id  SrSio8°a  5f°xr  '  !ii,r«n. seismographs  was  to  detect,  S 
30  sec  „u  ratio,  a  5  x  10  ^  strain  in  a  narrow  band  centered  on 

pressure  induced^o^^  “J[a'as  microseismic.n°ise  and  residual  barometric 

met  with  an  extra “of ^  been 

son  between  the  S325L1  and  S325L2  seismograL  ■  The  C°mpari- 

coincide  from  about  12  to  20  sec.  g  me  6325L1  and  S325L2  spectra 


This  observation  is  significant.  It  means: 


S32S-1  'strain  s^s^™; C°”Stant  “  Wn  ‘°  a=FUracy  on  the 


j  ±  1.96  dB 

90%  CONFIDENCE 
-1  LIMIT 

3841  SAMPLES 
2S6  LAGS 
2  SPS 

PARZEN  SMOOTHING 


2.  The  calibration  constant  is  known  to  sufficient  accuracy  on  the 
S325-2  strain  seismometer; 

3.  In  the  12  to  20  sec  period  range  both  strain  seismometers  are  seeing 
the  same  identical  ground  strain  which  must  be  a  true  representation  of  the 
ground  motion. 

Another  observation  to  be  made  from  these  spectra  is  that  the  S325L2  has  more 
power  than  the  S325L1  at  periods  longer  than  about  25  sec  and  has  less  power 
at  periods  shorter  than  about  12  sec  including  the  8  sec  microseismic  peak. 

This  same  frequency  dependent  phenomena  will  be  seen  in  the  other  comparisons 
that  follow.  The  true  mechanism  acting  here  is  not  positively  known.  However, 
a  working  hypothesis  is  presented  that 

a.  Both  the  S325-1  and  S325-2  are  attenuated  to  signals  from  travelling 
waves  in  the  e^rth  by  a  sma^l  factor  (see  section  11.4.4)  either  because  of 
the  joint  and  fault  system  in  the  rock  in  the  vicinity  of  the  mine  or  because 
of  the  geometry  of  Goldmine  Mountain  and  the  surrounding  region  and; 

b.  The  S325L2  is  acted  on  by  an  additional  frequency  dependent  damping 
force  proportional  to  the  rate-of-change  of  ground  strain.  The  S325-1  strain 
seismometer  is  mounted  in  solid  rock  both  at  the  fixed  end  anchor  and  at  the 
transducer,  while  the  S325-2  strain  seismometer  is  mounted  on  a  different 
mineralized  fault  zone  at  both  ends.  The  two  faults  both  dip  at  about  70  deg. 

The  fault  at  the  northern  end  dips  to  the  northeast  and  the  fault  at  the' 
southern  end  dips  to  the  southeast. 

It  is  possible  that  the  small  attenuation  factor  may  be  different  for  both 
instruments.  If  it  were  the  same,  the  two  spectra  in  figure  20  would  be  equal 
at  long  periods  and  the  S325L2  would  have  less  power  than  the  S325L1  at  all 
periods.  The  attenuation  factor  could  be  the  same,  and  the  small  difference 
at  the  longest  periods  in  figure  20  could  be  a  calibration  uncertainty. 

The  coherence  between  the  two  parallel  strain  seismometers  peaks  at  the  microseismic 
periods  of  8  sec  and  18  sec.  While  the  coherence  of  about  0.25  at  30  sec  and 
from  50  to  256  sec  is  not  statistically  significant,  it  is  felt  that  the 
coherence  of  a  partially  deterministic  signal  plus  noise  is  a  more  robust 
process  than  indicated  by  statistical  tests  with  random  numbers  and  that  this 
coherence  is  meaningful.  The  coherence  is  attributed  to  ground  strains  from 
the  small  to  moderate  microbarometric  pressures.  The  phase  angle  lead  of  the 
S325L2  over  the  S32Stl  is  seen  to  be  almost  zero  between  8  and  30  sec  and  30 
and  50  sec  and  almost  90  deg  lead  at  30  sec  and  between  50  and  256  sec.  The  , 

90  deg  phase  lead  is  in  agreement  with  the  above  hypothesis  that  the  S325L2 
is  acted  upon  by  a  damping  force  proportional  to  r£.te-of-change  of  strain. 

In  studying  the  numerical  values  of  the  two  power  spectral  densities  small 
peaks  not  discernible  in  the  plot  in  figure  21  were  noticed  at  adjacent  fre¬ 
quencies  for  the  S 325 LI  and  S325L2.  If  these  peaks  are  reaT7  the  analysis 
plotted  in  figure  21  is  a  true  representation.  However,  since  the  variations 
in  the  magnitudes  of  the  smaller  peaks  are  much  less  than  the  statistical 
stability  of  the  analysis,  they  may  not  be  real  physical  phenomena.  To  test 
the  validity  of  their  frequency  separation,  another  analysis  was  performed  with 
half  as  many  lags  (128  instead  of  256)  to  obtain  a  smaller  confidence  limit. 

The  results  of  this  analysis  are  plotted  in  figure  22.  These  results  are  seen 
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to  be  a  smoothed  version  of  figure  21  with  these  primary  differences.  The 
cross  spectral  peak  at  30  sec  in  figure  21  is  smoothed  into  the  peak  at  16 
sec  in  figure  22.  The  smoothing  improved  the  coherence  at  all  periods  in 
figure  22  over  figure  21.  The  phase  lead  of  S325L2  over  S325L1  is  almost  zero 
from  1  sec  to  about  25  sec,  then  smoothly  changes  to  a  90  deg  phase  lead. 

The  conclusions  arrived  at  from  the  analysis  of  figure  21  with  greater  frequency 
resolution  are  substantiated  by  the  analysis  of  figure  22  with  greater 
statistical  stability. 

/ 

The  computational  noise  referred  to  in  section  11.4.1  is  seen  to  be  on  the 
order  of  10_ 27  strain2/Hz  compared  to  maximum  values  on  the  order  of  10*22 
strain2/Hz  in  both  figures  21  and  22.  A  sharp  break  is  seen  in  the  slope  of 
all  of  the  spectra  at  about  6  sec.  At  periods  shorter  than  6  sec,  the  results 
should  be  considered  invalid. 

The  S325L1  strain  seismograph  spectra  are  compared  with  the  P325L  inertial 
seismograph  spectra  in  figure  23.  The  strain  spectra  and  the  strain-inertial 
cross  spectra  are  normalized  as  described  in  section  11.4.1  with  this  normali- 
>  zation,  the  two  sets  of  spectra  of  this  sample  are  almost  equal  at  25  sec. 

As  stated  above,  the  inertial  seismometer  had  been  recently  worked  on  and  was 
experiencing  long-period  noise  oscillations.  Consequently,  the  P325L  spectra 
at  periods  longer  than  about  40  Sec  should  be  used  only  for  qualitative 
comparisons.  Even  with  the  extra  long-period  noise,  the  total  mean  square 
displacement  (relative  to  25  sec)  in  the  passband  of  the  seismograph  was 
3.10  x  10*15  m2  which  is  55.6  mu  rms  in  the  horizontal  direction  of  325  deg 
azimuth.  The  three  peaks  in  the  cross  spectpa  are  significant.  The  two  peaks 
at  g  sec  and  16  sec  are  from  microseisms  and  the  10  sec  peak  is  from  some 
unidentified  signal.  The  significance  of  the  10  sec  peak  is  its  confirmation 
as  a  true  signal  even  though  it  was  not  apparent  in  the  spectra  of  the  S325L2 
seismograph  in  figure  21  or  22  and  is  almost  smoothed  out  of  the  P325L  spectra. 
The  significance  of  the  8  sec  and  16  sec  peaks  are  that  the  strain  and  inertial 
seismographs  do  not  differ  at  these  periods  by  the  normalization  factor 
applicable  for  the  calibration  period  of  25  sec.  The  fact  that  there  is  less 
energy  on  the  S325L1  strain  seismograph  than  the  P325L  inertial  seismograph 
implies  that: 

t  * 

a.  The  apparent  pjiase  velocity  of  this  energy  in  the  horizontal  plane 
in  the  325  deg  azimuth  is  larger  at  8  sec  and  16  sec  than  at  25  sec,  or 

b.  The  strain  in  the  325  deg  azimuth  is  not  coupled  to  travelling  waves 
in  the  Earth  to  ^he  same  degree  that  the  output  of  a  horizontal  pendulum  is 
coupled  to  the  same  waves,  or 

c.  The  local  geometry  of  Goldmine  Mountain  and  the  surrounding  rfegion 
has  modified  the  component  of  the  local  vector  wave  number  in  the  325  deg 
azimuth,  or 

d.  A  combination  of  these  factors . 

There  are  indications  that  all  of  these  have  contributed.  Signals  appear 
smaller  than  expected  on  the  two  parallel  325  deg  azimuth  seismographs  as  will 
be  shown  in  section  11.4.4.  There  also  is  a  definite  difference  in  the  apparent 
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phase  velocity  of  the  energy  at  8  sec  and  of  the  energy  at  16  sec  with  relatively 
less  16  sec  energy  being  seen  in  the  spectra  and  in  the  film  recording^. 

The  coherence  between  the  P325L  and  S325L1  plotted  in  figure  23b  suffered  from 
computational  noise  resulting  from  the  large  scale  difference  between  strain 
and  pendulum  data.  Because  of  the  noise,  the  data  should  be  interpreted 
qualitatively  as  high  coherence  at  8  sec  and  fairly  constant  coherence  at 
about  0.25  from  10  sec  to  256  sec.  The  phase  lead  of  P325L  is  seen  to  be 
zero  for  periods  Up  to  10  sec  and  close  to  180  deg  for  the  longer  periods. 

Figure  24  is  an  airalysis  of  the  same  data  used  for  figure  23  except  that  128 
lags  were  used  for  a  smaller  confidence  limit.  The  results  are  the  same  as 
discussed  above  except  th*t  the  10  sec  peak  was  smoothed  into  the  adjacent' 
passbands.  The  coherence  and  relative  phase  angles  are  the  same. 

r  c. ' 

The  perpendicular  strain  seismograph  S55L  is  compared  to  the  S325L1  in  figure  25. 
The  total  mean  square  strain  (relative  to  25  sec)  on  the  S55L  in  the  total 
passband  was  9.79  x  1C“25  strain^  which  is  an  rms  strain  of  9.89  x  10~13.  The 
two  spectra  are  very  similar  at  periods  up  to  about  20  sec  with  peaks  in  the 
cross  spectra  at  8,  10,  18  and  32  sec.  The  fifst  three  peaks  are  from  coherent 
energy.  The  8  sec  and  18  sec  peaks  and  most  periods  have  a  180  deg  phase 
relationship.  The  180  deg  phase  relationship  is  characteristic  of  transverse 
shear  strain.  At  the  10  sec  peak,  S55L  leads  S325L1  by  300  deg  and  at  32  sec 
by  240  deg.  Since  the  phase  angle  isi  not  180  deg,  the  10  and  32  sec  peaks  pro¬ 
bably  are  not  the  result;  of  shear  strains.  The  occurrence  of  the  30-32  sec 
peak  in  the  cross  spectra  of  both  the  parallel  strain  seismographs  and  the  per¬ 
pendicular  strain  seismographs  indicates  that  there  was  real  ground  motion  at 
this  period.  This  may  be  another  frequency  of  low-level  microseismic  activity. 
The  level  of  the  nower  sDectra  in  both  azimuths  is  relatively  the  same  with  the 

«  ii  * 

differences  well  within  the  range  of  possible  variation  in  the  azimuthal  response 
of  the  strain  systems  to  the  sources  of  microseismic  noise. 

O  „ 

11.4.3  High  Wind  Conditions 


During  episodes  of  high  wind,  there  are  variations  in  the  pressure  field  which 
create  displacements,  tilts,  and  strains  within  the  ground.  These  disturbances 

«re  discussed  in  more  detail  in  sections  11.5  and  11.6.  The  two  parallel 
orizontal  strain  seismographs  S325L1  and  S325L2  respond  identically  to  these  ' 
ground  strains.  Figure  26  is  a  portion  of  the  time  sample  that  was  analyzed 
for  a  comparison  during  windy  conditions.  The  wind  was  blowing  about  20  mph 
and  gusting  to  35  mph.  The  ML2  microbarograph  outside  the  mine  was  responding 
to  large  pressure  fluctuations.  In  the  figure  the  equal  and  opposite  strain 
typical  of  a  shear  stress  is  visible  at  several  points,  with  an  especially 
large  yet  typical  pulse  at  about  0503Z.  The  S325L1  and  S325L2  seismographs 
track  each  other  identically  in  the  time  domain  in  figure  26  and  the  spectra  * 
almost  overlay  in  the  frequency  domain  in  figure  27.  During  this  sample  the 
total  mean  square  strain  (referenced  to  25  sec)  in  the  passband  was  4.198  x 
10"23  strain^  for  the  S325L1  and  5.458  x  10-23  strain^  for  the  S325L2  and  the 
rms  values  were  6.479  10" 12  strain  ind  7.388  x  10_12  strain,  respectively. 

In  figure  27  the  spectra  are  seen  to  be  almost  identical  except  for  a  slight 
increase  in  value  for  the  S325L2  at  increasing  period.  The  coherence  is  one 
and  there  is  no  relative  phase  shift  between  the  seismographs.  The  seismo¬ 
graphs  are  recording  identically.  The  time  sample  was  extended  to  6.7  fir) 
and  analyzed  for  the  comparison  for  periods  between  20  and  2560  sec  seen  in 
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Power  spectral  density,  coherence,  and  relative  phase  angle  wi 

P325LaKer?iarrH  1  h?rizontdl  seismographs,  S325L1  strain  a 
P325L  inertial,  during  low  wind  conditions  c 


TR  70-29 


PHASE  (dagrwn)  COHERENCE  POWER  SPECTRAL  DENSITY  <<STRAIN)2/Hz  RELATIVE  TO  25  SEC,  SEE  TEXT) 


Power  spectral 
density 


PERIOD  (neon*) 


Coherence 


Phase  lead 
S55L  over  S325L1 


PERIOD  (tacondt) 


Figure  25.  Power  spectral  density,  coherence,  and  relative  phase  angle  with 
256  lags  for  perpendicular  horizontal  strain  seismographs  S325L1 
and  S55L  during  low  wind  conditions 
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figure  28.  The  results  are  the  sane  as  at  the  shorter  periods  except  between 
300  and  1280  sec  where  the  S325L1  has  more  power  than  the  S525L2.  The 
coherence  is  almost  one  except  in  the  anomalous  period  range  and  the  relative 
phase  angle  remains  at  essentially  zero  for  all  periods.  A  comparison  between 
the  ULP  response  seismographs  S325U1  and  S325U2  is  illustrated  in  figure  29. 

The  rms  strains  (referred  to  25  sec)  during  the  6.7  hr  period  in  the  total 
passband  of  these  seismographs  were  3.61  x  10~H  strain  and  2.08  x  10-H  strain, 
respectively.  The  two  a  ?  seen  to  be  essentially  the  same  at  periods  out  to 
234  sec  and  then  the  S32.>U1  exhibits  considerably  more  power  than  the  S325U2. 

During  the  pressure  induced  disturbances,  the  horizontal  inertial  seismometers 
are  subjected  to  considerable  tilting  as  can  be  seen  in  the  seismogram  in 
figure  26  on  the  inertial  ALPS  response  channels  (P-L).  The  tilts  are  more 
evident  on  the  inertial  extended-long-period  response  (P-X)  which  has  a  higher 
magnification  at  the  longer  periods  than  the  ALPS  response.  The  inertial  P525L 
is  compared  with  the  strain  S325L1  in  figure  30.  The  rms  ground  displacement 
(referred  to  25  sec)  during  this  32.0  min  sample  was  33.0  mu.  This  value  is 
lower  than  the  rms  displacement  during  the  low  wind  sample  because  in  the 
intervening  37  hr  the  seismometer  settled  down  to  almost  normal  operation. 

In  the  P325L  spectra  in  figure  30,  the  general  background  spectrum  hds  risen 
until  the  microseismic  peaks  at  8.4  and  16.2  sec  are  almost  swamped  by  the 
general  rise  in  the  spectrum.  The  coherence  is  relatively  low  indicating  that 
the  tilt  response  of  the  horizontal  inertial  seismometer  is  only  partially 
linearly  related  to  the  ground  strain.  The  ground  strain  has  increased  so 
that  the  entire  spectra  is  above  the  strain  of  the  microseismic  peaks  which 
are  not  visible  either  in  the  analog  trace  or  as  a  spectral  peak. 

The  perpendicular  horizontal  strain  seismographs  are  very  comparable  in  spectra. 
The  equal  amplitude  and  180  deg  phase  relationship  is  indicative  of  the  res¬ 
ponse  of  perpendicular  linear  strainmeters  to  transverse  shear  strain.  During 
the  high  wind  sample  analyzed,  the  rms  strain  (referenced  to  25  sec)  for  the 
S55L  seismograph  was  6.500  x  10‘12  strain.  This  rms  value  is  within  0.32 
percent  of  the  6.479  x  10~12  rms  strain  on  the  S325L1.  The  fact  that  the  rms 
values  are  this  close  is  fortuitous  considering  the  errors  in  calibration  and 
analysis,  but  it  does  substantiate  that  calibrations  are  very  close  and  are 
well  within  acceptable  operational  tolerances.  The  spectral  comparison  is 
plotted  in  figure  31.  The  spectra  at  periods  shorter  than  about  10  sec  are  in 
the  computational  noise  and  are  not  meaningful.  The  spectra  are  seen  to  be 
very  smooth  and  devoid  of  any  peaks.  The  cross  spectra  does  show  a  minor  peak 
at  the  microseismic  period  of  16  sec.  This  small  amount  of  microseismic 
strain  apparently  mixes  with  the  pressure  induced  strain  in  a  nonlinear  manner 
resulting  in  the  dip  in  the  spectra  at  20  sec.  The  coherence  holds  to  reasonable 
values  except  for  a  minimum  at  20  sec  and  at  periods  longer  than  100  sec.  The 
phase  angle  is  almost  180  deg  at  all  periods.  Figure  32  extends  the  comparison 
to  longer  periods  with  the  ULP  response  seismographs  S325U1  and  S55U.  About 
the  same  picture  is  seen  as  before.  The  spectra  are  similar  in  appearance  out 
to  periods  near  250  sec.  A  comparison  of  figures  29  and  32  shows  that  at  the 
longer  periods,  the  S55U  resembles  the  S325U2  more  than  the  S325U1.  The  cross 
spectra  have  many  sharp  high  Q  peaks.  With  the  frequency  resolution  and 
statistical  stability  of  the  analysis,  these  peaks  are  considered  as  being 
real  and  are  definitely  not  an  artifact  of  the  data  analysi£-_^Also ,  since  these 
peaks  appear  in  the  cross  spectra  of  the  perpendicular  s*rains/\but  not  in  the 
cross  spectra  of  the  parallel  strains,  they  must  be  associated  w\th  a  shear  strain 
phenomena. 
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Figure  28. 


Power  spectral  density,  coherence,  and  relative  phase  angle  of 
parallel  horizontal  strain  seismographs  S325L1  and  S325L2  during 
high  wind  conditions,  20  to  2560  sec  range  K 
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Figure  29.  Power  spectral  density,  coherence,  and  relative  phase  angle  of 

parallel  horizontal  strain  seismographs  S325U1  and  S325U2  during 
high  wind  conditions 
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Figure  30.  Power  spectral  density,  coherence,  and  rela 
of  parallel  horizontal  seismographs,  S325L1 
P325L  inertial,  during  high  wind  conditions 
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Figure  31.  Power  spectral  density,  coherence,  and  relative  phase  angle  of 
perpendicular  horizontal  strain  seismographs  S325L1  and  S55L 
during  high  wind  conditions 
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coherence  of  the  peaks  in  the  cross  s  t  ,SeismoSraPh<5  shows  the 
relative  phase  angle  tends  toward  180^  between  about  0.25  and  0.4  The 

and  tends  toward  0  deg  at  oth"r  penodf  8  Peri°ds  be,“e"  and  210  Uc 


11.4.4 


Earthquake  Signals 


seismographs  producfdat1he''s^eUoutputs Vfo^t  ^  parallel  horizontal  strain 
are  associated  with  the  large  Peru  earth°  l raV^lllng  wav«s .  All  three  signal 

coLat0oftMS  ^3rthquake  are:  0  =  2023:2?Ufe9°2s3178a8IV19h0‘  USC5GS  PDE  Card 

t  of  northern  Peru,  =  6.6,  Ms  =  7  8  ’&  -%?  <r  j8'*  h  “  43  km>  near 
epicenter  137.4  deg.  The  first  /'8’  A  '  52 • 5  -deg,  azimuth  station-to- 

from  the  coda  of  the  LR1  Rayleigh  wave  during  75*  ^  f5  ebrthclua^  was 
group  velocities  of  1.72  to  1  10  km/ser  n,8  tlme  ulndow  equivalent  to 

S'gnal  fro.  the  LR2  Rayleigh  wtrvc  with  a  grau/'ew,  Slgnal  Saniple  *S  a  “"bi"ed 
2.20  km/sec  and  from  the  103  invP  j  ^p  velocity  wmtow  of  2.52  to 

velocity  window  of  3.37  tJ  2.95  kJSec  n  7  RafIeigh  “a«  "ith  .group 
of  high  gam  channels  were  driven  off  the  fill!8,  j,h  °f  theSe  time  sa«ples  all 
tape  recorders.  Therefore,  the  low  gain  and  were  clipping  the  magnetic 

used  for  the  comparisons.  The  third  earthni  Seismo8raPh  channels  were 
the  LR1  Rayleigh  wave  from  an  aftershock  ^e6^18"3*  1S  3  ma-’°r  Portion  of 

betWe6n  3'37  and  2‘95  Wsec  UsSgs'V^F  C™SPonds  to  group 

aftershock  are:  0  =  0136:10.2  9  3S  7Q  nw  u  _ ?GS,  PDE  Card  data  for  the 
northern  Peru,  n^,  =  6.0,  M<.  =  5.5'  A’_  r\  \  4?  km>  off  the  coast  of 

138  deg.  The  high  gain  seismographs  were  used’ ^  station-t°-epicenter  = 
sample  even  though  the  S325L1  andPS325f 7 J *  f  ^  analysis  of  the  third 
0203Z  and  02062  and  the  P325L  was  clin  pH  w  lpped  on  tbe  Peaks  between 
Tbe  total  time  of  clipping  was  about  ?P  d  ^  the  peaKs  from  0203Z  to  0207Z 

did  "0^  materially  S'Tat  th^  7°le  ^  W1"d°"'  a"d 

spectrum.  y  except  at  the  short-period  end  of  the 

r  f:re  33- 

channels  stopped  clipping.  Thus  this  sairroleV^^  after  the  low  gain  tape 
that  can  be  recorded  without  distortion  on  thelupl  S  ^  maximuni  signal 
are  recorded  at  two  signal  levels  J£"  7  ™  S  resP0nse  channels  which 
with  either  a  little  less  po^  ™  th^  sSuTat^  ^  be  aI“>st  identical 
than  20  sec  or  a  slight  calibration  errof  fte  Lhe  “  M  periods 
relative  phase  angle  is  zero  at  all  a  coherence  ls  one  and  the 

O"  J°  25  ^c),  in  the  whole  p  Jsba^fwas^"'?!'8* 'Jo-lO3"’’16, the  ™S  strai" 
S325LU  and  1.91  x  10-10  strain  on  the  S3251L2.  strain  °"  the 

during  the  portion'of "he'second^smnleTsh6*"  -°  p'’3'1"  each  other  identically 
figure  35,  show  the  same  results  fo7the  S375.  m  fH,34'  ^  sP«trd, 

S32oLLl  either  has  slightly  less  nower  n  ELI  and  S325LL2.  Again  the 

coherence  is  one  at  all  periods  except  10  sec”3^  ^anbration  error.  The 
relative  phase  angle  is  zero.  The  spectra  of  ^  ^  ^  n°iSe  ^  the 

noise  is  evident  as  the  approximately  whifp  f  k  5  pG  recorder- reproducer 
analog  antialias  filter  acting  uport  the  tanp"01-6  between  3  and  10  sec.  The 
sec.  The  same  antialias  filtfr  wai  used ?£' ",n01Se  15  clear  between  1  and  3 

spectra  were  calculated  in  the  1  to  256  sec  pafs'banT8  ^  th®  d3ta  f°r  Whicb 

sec  passband.  The  action  of  the  filter 
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Figure  33.  Power  spectral  density,  coherence,  and  relative  phase  angle  of 
parallel  strain  seismographs  S325LL1  and  S325LL2  during  the 
LR1  Rayleigh  wave  coda  of  the  Peru  earthquake 
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Figure  35.  Power  spectral  density,  coherence,  and  relative  phase  angle  of 
parallel  horizontal  strain  seismographs  S32SLL1  and  S325LL2 
during  the  LQ3  Love  wave  and  the  LR2  and  the  LR3  Rayleigh  waves 
of  the  Peru  earthquake 
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is  not  seen  in  any  of  the  other  figures  because  the  spectra  between  1  and  5 
sec  in  all  the  other  figures  is  well  below  the  computational  noise  level. 

The  peak  of  the  spectra  at  about  10*18  strain^/Hz  in  figure  33  and  the  tape 
noise  at  about  10*24  strain2/Hz  in  figure  35  indicate  a  potential  60  dB  dynamic 
range  with  narrow  band-pass  spectral  analysis.  The  same  dB  value  applies  to 
both  power  and  amplitude. 

The  strain  spectra  from  the  perpendicular  S325LL1  and  S55LL  are  plotted  in 
figure  36.  The  spectra  are  similarly  shaped  but  the  S325LL1  has  only  about 
two-thirds  the  magnitude  of  the  power  on  the  S55LL.  This  difference  is 
attributed  to  a  possible  cancellation  of  LR2  and  LR3  Rayleigh  motion  on  the 
almost  radial  325  deg  azimuth  instruments.  The  coherence  is  good  and  the 
relative  phase  angle  is  180  deg.  The  rms  strain  (referred  to  25  sec)  in  the 
whole  passband  for  this  sample  was:  S325LL1  1.500  x  10"H  strain,  S325LL2 

1.704  x  10*11  strain,  and  S55LL  1.983  x  10"H  strain. 

The  third  earthquake  signal  further  confirms  that  the  two  parallel  strain 
seismographs  give  the  same  output.  The  analog  signals  on  the  low  gain  channels 
are  seen  tc  track  identically  in  figute  37.  The  spectra  essentially  overlay 
in  figure  38.  The  data  below  7  sec  are  into  the  noise.  The  peaks  at  the 
shorter  periods  are  the  result  of  using  clipped  data  for  about  5  percent  of  the 
sample.  The  coherence  is  one  and  the  relative  phase  angle  is  zero  at  all 
periods.  Because  of  the  clipped  data  the  spectra  were  also  calculated  with 
128  lags  for  better  stability.  The  results  plotted  in  figure  39  are  the  same 
as  in  figure  38  except  that  the  noise  only  affects  periods  shorter  than  4  sec. 
The  coherence  is  one  at  all  periods  indicating  a  more  reliable  statistical 
estimate. 

Figure  40  contains  the  spectra  of  the  S325L1  and  the  P325L  seismographs.  There 
are  two  significant  differences.  First,  at  25  sec,  the  S325L1  is  a  factor  of 
3.5  lower  in  (equivalent)  power  or  1.9  lower  in  amplitude  than  the  P525L. 
Second,  the  cross  spectra  very  closely  follows  th?  S325L1  in  shape  and  power. 
There  are  two  interpretations  that  can  be  taken  from  these  data.  First,  the 
similarity  of  the  strain  spectra  and  the  cross  spectra  suggest  that  the  strain 
seismograph  S325L1  represents  the  linear  portion  of  true  ground  motion  and  the 
additional  power  on  the  inertial  seismograph  P325L  represents  a  nonlinear 
output  generated  by  localized  tilting  of  the  nearby  rock.  The  concept  of 
nonlinear  tilt  response  of  the  rock  is  discussed  in  section  11.6. 

The  second  interpretation  is  that  the  325  deg  azimuth  component  of  the  vector 
wave  number  in  the  vicinity  of  the  instruments  has  been  altered  from  the 
expected  wave  number  component.  Therefore,  a  better  comparsion  could  be  made 
by  emperically  adjusting  the  strain  seismograph  output  by  the  3.5  power  ratio. 
Even  when  this  adjustment  is  made,  the  inertial  seismograph  has  a  significantly 
greater  power  than  the  strain  seismograph  at  periods  longer  than  25  sec. 

A  nonlinear  response  to  tilting  could  smooth  out  the  spectral  peaks  on  the 
inertial  seismogram  also.  This  lowering  of  the  Q  could  result  from  a 
mechanical  friction  between  blocks  along  the  local  fault  and  joint  pattern. 

This  same  phenomena  may  be  causing  the  frequency  dependent  relationship 
between  the  two  parallel  strain  seismographs.  The  9  sec  peak  in  figure  40  is 
not  as  high  a  Q  on  the  P325L  as  it  is  on  the  S325L1. 
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Figure  36.  Power  spectral  density,  coherence,  and  relative  phase  angle  of 
perpendicular  horizontal  strain  seismographs  S325LL1  and  S55LL 
during  the  LQ3  Love  wave  and  the  LR2  and  the  LR3  Rayleigh  waves 
of  the  Peru  earthquake 
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Figure  38.  Power  spectral  density  coherence,  and  relative  phase  angle  with 
256  lags  for  parallel  horizontal  strain  seismographs  S325L1  and 
S325L2  during  the  LR1  Rayleigh  wave  of  an  aftershock  of  the 
Peru  earthquake  G6(J 
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Figure  40.  PowOr  spectral  density,  coherence,  and  relative  phase  angle  with 
128  lags  for  parallel  horizontal  seismographs,  S325L1  strain  and 
P325L  inertial,  during  the  LR1  Rayleigh  wave  of  an  aftershock  of 
the  Peru  earthquake.  G6021 
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In  figure  40,  the  coherence  is  about  0  9  near  8  see  _ 

phase  angle  is  zero  except  at  periods  of  64  d  “u  SeC'  reldtlve 

During  the  third  earthquake  signal  anaU^ed  rt  *eC  11  is  180  de8‘ 

sec)  in  the  whole  passband  was-  S325Ll'  1  82  /i  n*  '  i"1011011  (relative  to  25 
10"11  strain,  and  ?325L  9.05  x  10  7  J  ' 82  X  10  Strain’  S325L2  ^  x 

11-5  SEISMIC  NOISE  RESULTING  FROM  ATMOSPHERIC  PRESSURE  VARIATIONS 
been^knoim^for'severa^years^ftecent  1  ?°i*)r  "lows  has 

created  by  wiati^s  ij  the  ?P  resuUl"8  fr<»  forces 

good  success.  This  section  of  rheP  ^ C  Presfure  field  can  be  predicted  with 

of  pressure  induced  noise-  extends  the^ec^t l"")  °f  the  past  observations 

ground  displacements  to  include  ground  strain ^„d”i  t^and’''1”’  l°  prediCt 
comparison  between  the  strain  cPkmn  .  ,  ,  tl?t’  and  presents  a  spectral 

11-6  extends  the  noise  considerations^  include  Hneartnd°n^P?:  SeCti°n 
on  horizontal  inertial  seismoeranhs  a ■  7  r  and  non-linear  tilt  noise 

5Hrr ■"--"-noise 

Some  of*  the  motions  reft? rrpH  tr>  in  t-u *  , 

Anpstrnmc  wu •  1  -  *  in  two  sections  are  on  the  order  1  to  n  i 

i^:T^besin  to  detect  the 


11.5.1 


Previous_Observations  of  Wind-Related  Long-Period  Noise 


dandcS:ednhet™K";  ias  ta  1 1  at  ion  eto‘ in  crease  t'he'i  rTtZufy  will'™;""'* 
summary  of  long-peri.™ I^gr^ ^ 

a^SSSSSs^ 
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Under  the  same  author's  (JEF)  supervision,  Milam  (1965)  continued  the  development 
of  improved  installation  techniques  to  minimize  all  noise  in  the  ALPS 
system.  Pressure  induced  tilt  was  reduced  by  a  fuither  significant  amount 
on  a  horizontal  seismometer  installed  in  an  8  m  deep  vault  about  100  m  from 
the  walk-in  valut  and  the  1  m  deep  tank  vault.  A  secord  seismograph  installation 
was  made  about  100  m  from  the  surface  in  the  Stevensor -Bennett  Mine  near  Las 
Cruces,  New  Mexico  (LC-NM).  (A  Long-Range  Seismic  Measurements  (LRSM)  Program 
team  co-occupied  the  mine  during  these  experiments  but  operated  independentl) .) 
The  Las  Cruces  installation  verified  the  capability  of  the  full  three -component 
seismograph  to  operate  at  magnifications  of  100K  (at  X10  view)  for  extended 
periods  of  time.  During  this  experiment,  the  need  to  use  tank  vaults  inside 
a  sealed  seismometer  chamber  well  within  a  mine  was  conclusively  established. 
After  all  pressure  and  temperature  effects  acting  diTectly  on  the  seismometers 
were  eliminated,  noise  was  observed  on  the  horizontal  seismometers  when  the 
wind  velocity  exceeded  about  30  mph .  The  existence  of  a  wind  related  noise 
threshold  at  both  KM0  and  LC-NM  confirms  that  the  elastic  loading  is  pro¬ 
portional  to  the  pressure  induced  force  at  the  surface.  Since  the  wind  velocity 
related  to  the  noise  threshold  at  a  100  m  depth  at  Las  Cruces  was  greater 
than  the  wind  velocity  related  to  the  noise  threshold  in  the  1  m  deep  or  the 
8  m  deep  tank  vaults  at  WM0,  the  observations  confirmed  that  the  amplitude 
of  the  displacements  was  related  to  the  wind  velocity. 

A  s  milar  attenuation  of  pressure  induced  noise  was  observed  at  the  l\ MO  short- 
period  strain  installation  in  1967.  (See  Shopland,  1968.)  The  original 
trench  depths  were  2.8  to  3.7  m  on  the  east  strain  seismometer  and  3.6  to 
3.7  m  on  the  north  strain  seismometer.  An  additional  1.3  m  of  overburden 
was  added.  The  wind  velocity  related  to  the  noise  threshold  increased  from 
about  7  km/hr  on  the  east  seismograph  and  about  20  km/hr  on  the  north  seismo¬ 
graph  to  about  35  km/hr  on  both  seismographs.  A  spectral  comparison  of  the 
east  seismograph  during  15  to  20  mph  winds  before  and  after  the  additional 
overburden  was  added  showed  a  5  to  ,10  dB  power  reduction  at  frequencies  above 
0.8  Hz  and  not  much  change  in  power  at  longer  periods.  Figure  41  shows  a 
replot  of  the  data  points  in  figure  51  from  Shopland  (1968)  on  a  log- log  scale. 
Straight  lines  have  been  visually  fit  to  these  data  points.  The  equations 
for  the  lines  are: 

strain  east  before 

u  =  0.1625  C1'495 


strain  north  before 


u  =  0.000324  C2-91 

strain  north  after 


u  =  0.000563  C2-20 


where 


u  =  equivalent  ground  motion  in  millimicrons 
peak-to-peak  calculated  from  the  strain 
seismograms 

C  =  wind  velocity  in  km/hr 
,  -68- 
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Figure  41. 


Relationship  between  equivalent  ground  displacement  and 
wind  velocity  at  WMO 
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The  improvement  of  the  operation  of  a  long-period  seismograph  by  burial  in  a 
culturally  noisy  location  was  also  demonstrated  by  Rector  (1965).  Three 
seismometers  in  a  15.25  m  deep  cased  hole,  in  a  surface  tank  vault  inside  a 
building,  and  on  a  concrete  pier  in  a  building  were  used  for  the  tests  at  the 
Geotech  Garland  facility.  After  properly  pressure  sealing  and  insulating 
the  buried  seismometer,  operating  magnifications  of  10  to  20K  were  possible. 
During  the  same  time  intervals,  the  two  surface  seismometers  could  only 
operate  at  0.5  to  2K.  The  burial  sufficiently  attenuated  most  cultural 
vibrational  and  pressure  induced  noises;  however,  all  seismometers  responded 
to  large  tilts  resulting  from  surface  loading  by  trains  and  vehicles.  No 
quantitative  tilt  comparisons  could  be  made  because  of  the  relative  locations 
and  depths  of  the  various  piers. 

Other  cases  of  improved  long-period  seismograph  operation  were  reported  in 
the  LRSM  program. barren  (1966),  excavated  a  9.1  m  deep  vault  at  La  Paz, 

Bolivia.  In  conjunction  with  improvements  made  in  the  thermal  and  pressure 
stability  the  depth  of  the  installation  attenuated  kind  related  disturbances 
and  resulted  in  an  improvement  of  operating  magnifications  from  the  10  to 
20K  range  to  the  100  to  200K  range.  Simons  (1968)  reported  on  an  experiment 
that  compared  the  operation  of  long-period  seismographs  in  surface  tanks  and 
at  the  bottom  of  162  ft  deep  missile  silos  at  two  different  locations  both  .  ^ 
of  which  were  in  sediments.  He  concluded  that  there  was  noise  on  the  horizontal 
seismometers  that  positively  correlated  with  wind  gusts  (more  than  a  steady 
wind).  However,  no  quantitative  relationship  could  be  established  between 
tilt  at  the  surface  and  at  depth,  in  fact,  at  times  the  surface  horizontal^ 
seismometers  would  tilt  first  and  at  other  times  the  silo  seismometers  would 
tilt  first.-  The  relative  amplitudes  of  the  tilts  also  varied  with  either 
the  surface  or  the  silo  seismometers  having  the  largest  amplitude. 

The  Uinta  Basin  Observatory  (UBO)  operated  three  sets  of  three -component 
long-period  seismographs  simultaneously  (Teledyne  Industries,  1966).  Ine 
approximate  depths  to  the  three  piers  were  0.6  m,  2  m;  and  15  m.  By  visual 
comparison  of  recordings  made,  the  deeper  vaults  during  windy  conditions  were 
about  6  dB  and  about  18  dB  quieter,  respectively,  than  the  shallowest  vault. 

A  detailed  quantitative  study  of  the  relationship  between  the  wird  velocity, 
pressure,  and  seismic  noise  versus  depth  was  not  made. 

Z'"'  '  > 

During  testing  of  the  triaXial  long-period  seismometer  in  a  60  m  (200  ft)  deip 
hole  at  UBO,  the  coordinate  transformed  hofizcmtal  outputs  of  the  triaxia/. 
system  were  compared  to  the  three- component  system  at  the  15  m  depth.  The 
triaxial  system  output  from  a  depth  of  53  tn  (175.  ft)  was  found  to  be  abou 
12  dB  quieter  than  the  output  of  the  shallower  system  (Kirkpatrick,  1968.) 

H.5.2  Long-Period  Earth  Strains  and  Atmospheric  Pressu?>i^riatie>£ 

To  explain  the  many  observations  of  (1)  wind  related  noise  on 
seismographs  and  (2)  attenuation  of  this  noise  with  depth,  Sorrells  (1969)  , 

Der  (1969),  and  Sorrells  and  Der  (1970)  developed  a  theory  for  the  displace¬ 
ments  associated  with  the  passage  of  a  plane  wave  of  air  pressure  over  an 
infinite,  isotropic,  homogeneous,  half-space  and  over  a  layered  half-space . 
Their  equations  for  displacements  ip  the  half-space  have  been  differentiated 
to  obtain,  the  relationship  between  strains  and  tifts  and  the  amplitude  cf  th  r 
pressure  wave.  '  ?• 
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Consider  a  system  of  axes  in  which  the  1  axis  is  in  the  direction  of  wave 
propagation  and  the  3  axis  is  vertical  down.  The  displacements  derived  by 
Sorrells  (Sorrells  and  Der,  1970)  with  a  slight  change  in  nomenclature  are 


By  taking  the  appropriate  derivatives,  the  linear  strain  in  the  direction  of 
propagation,  the  vertical  linear  strain,  the  shear  strain  in  the  1-3  plime, 
and  the  tilt  in  the  1  direction  are  found  to  be: 


f 


! 


.  a‘  linear  strains  and  the  tilt  do  not  have  the  1/k  =  c„/l  '  in 

itr k  is  wave  7,115  »“»“■»  ^ 


mn 


(8) 


m  =  n . 


(9) 


Since  the  velocity  of  the  pressure  wave  is  much  less  than  the  velocity  of 
seisnuc  serface  waves,  k  is  nuch  larger  for  the  pressure  wave.  Consequently 
the  ratio  of  pressure  induced  noise  to  surface  wave  signal  on  the  strain  Y’ 

inertia^seismograph i0"a,,Y  8reater  ““  the  rati°  °f  "oi«  «  -8"al  on  the 


„  ,  b-  Tbe  linear  strains  are  either  in-phase  or  180  deg  out-of-phase  with 
each  other,  depending  on  the  sign  of  the  terms  in  the  parenthesis. 


e  isC;f,/°r  e3Ch  W3Ve  nUmber*  there  is  0ne  ^Pth  where  the  horizontal  strain 
^11  zero. 


periods  rh^  ^  phaSe  Velocity>  at  the  longer 

periods,  the  term  containing  x3  m  the  parenthesis  becomes  negligible. 


e-  1116  linear  strains  ejj  and  e33  are  ±90  deg  out  of  phase  with  the  n, 
displacement,  the  shear  strain  ei3,  and  the  tilt.  1 


and  infor  lfin  t**  °Ut  phaSe  with  the  Vertical  displacement  u3 

and  m  or  180  deg  out  of  phase  with  the  horizontal  displacement.  3 


g\  7,16  faCt°r  controlling  the  attenuation  with  depth  of  all  displace¬ 
ments,  strains,  and  tilt  is  the  exponential  term  e*kox3. 


The  amplitude  of  the  pressure  induced  noise  on  strain  and  inertial  instruments 
SrT  rapidly  with  depth  than  the  amplitude  of  surface  waves  decays 
with  depth  Therefore,  by  increasing  the  depth  of  an  installation,  the 
p  essure  induced  noise  can  be  reduced  for  better  performance.  This  improve- 
ment  in  signal -to-noise  ratio  is  one  of  the  main  reasons  why  a  mine  was 
elected  for  the  installation  site  on  this  project. 


11.5.3  Observed  Strains  and  Pressure  Variati 


ions 


^  secnln  u  gramS  r!COrded  durir*  the  hi®h  *ind  co^di^ions  discussed 

?  11-4-31were/omP^red  Wlth  ^e  ML2  microbarograph  outside  the  mine. 
The  result  are  plotted  m  figures  42  and  43  The  rms  values  (referred  to 

forSsLm  ?  naf^in-ll1^08^5  in  the  whole  Passband  were  3.605  x  1Q-H 
for  S325U1 ,  2.082  x  10  H  for  S325U2,  and  1.647  x  10'H  for  S55U.  .  The  rms 

va  ue  (referred  to  112  sec)  for  the  ML 2  microbarograph  was  22.02  ubar.  If  the 

total  strain  enepgy  on  all  three  seismographs  is  assumed  to  be  caused  by  the 

pressure  source,  then  the  ratio  of  rms  values  will  give  a  ratio  of  strain  per 
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microbar  disturbance.  Taking  these  ratios  gives  1.64  x  10'1-  strain/„bar, 
9.46  x  10*1^  strain/iibar,  and  7.48  x  10*12  strain/ubar,  respectively.  This 
assumption  is  not  correct  as  indicated  by  the  coherence  in  figures  42  and  45. 
The  spectra  in  both  figures  were  plotted  with  scales  arbitrarily  equalized 
at  1  ubar  =  10*12  strain.  A  more  meaningful  relationship  will  be  obtained  in 
future  work  by  dividing  the  cross  spectra  by  the  microbarograph  spectra.  The 
cross  spectra  is  an  estimate  of  the  linearly  related  power  in  the  microbaro¬ 
graph  and  the  seismograph  and  does  not  include  microseismic  power  or  system 
noise  power  from  either  instrument. 


11.6  LINEAR  AND  NON-LINEAR  TILT  NOISE  ON  LONG-PERIOD  SEISMOGRAPHS 

The  ground  tilts  with  the  passage  of  Rayleigh  waves  and  with  the  application 
of  localized  forces.  Seismographs  respond  to  these  tilts  with  varying  degrees 
of  sensitivity.  Linear  strain  seismographs  respond  to  differential  motion 
along  their  sensitive  axis  and  are  insensitive  to  motions  orthogonal  to  this 
axis.  Therefore,  a  linear  strain  seismograph  constructed  to  measure  the 
strain  component  e^  will  be  insensitive  to  the  other  eight  strain  components 
ejj ,  i  =  1,2,3,  j  =  1,2,3,  except  i  =  j  =  1.  The  vertical  inertial  seismograph 
will  respond  to  the  vertical  displacement  component  of  the  tilt.  Generally 
for  pressure  induced  noise,  the  output  of  a  vertical  seismograph  is  about 
1/10  the  output  of  the  horizontal  seismograph.  The  horizontal  seismograph  must 
respond  in  the  same  manner  to  a  true  horizontal  acceleration  as  it  does  to 
the  component  along  the  sensitive  axis  of  the  instrument  of  the  gravitational 
acceleration  created  by  a  tilt  of  the  base  on  which  the  seismometer  rests. 

This  part  of  the  report  will  (1)  adapt  a  general  theory  for  non-sensitive 
axis  di  ‘jrbances  to  horizontal  pendulum  seismometers  developed  by  Rodgers 
(1968)  to  the  specific  parameters  of  the  ALPS  response  for  tilt  resulting 
from  Rayleigh  waves  and  pressure  loading,  (2)  present  some  data  on  the  nature 
of  the  change  of  the  seismic  spectrum  with  increasing  wind  velocity  and  pressure 
variations,  and  (3)  present  data  that  indicate  much  of  the  tilting  observed 
in  horizontal  inertial  seismographs  is  the  result  of  rigid  block  non-linear 
tilting  and  amplification  of  the  tilts  significantly  above  predicted  values. 

11.6.1  Tilt  Theory  for  ALPS  Response  Horizontal  Seismographs 

Rodgers  (1968)  very  thoroughly  treats  the  response  of  the  horizontal  pendulum 
seismometer  to  Rayleigh  and  Love  waves,  tilts,  and  free  oscillations.  The 
most  significant  spurious  response  is  to  tilt  and  to  the  tilt  associated  with 
Rayleigh  waves.  These  two  spurious  responses  will  be  considered  in  this  para¬ 
graph.  As  mentioned  above,  the  horizontal  pendulum  responds  to  the  component 
of  the  gravity  acceleration  created  by  the  tilt  of  the  base.  Thus  the  seis¬ 
mometer  responds  to  tilt  as  an  accelerometer.  If  Ap  (u)  is  the  displacement 
amplitude  frequency  response  and  Aj  (u)  is  the  tilt  amplitude  frequency,  with 
g  the  gravitational  acceleration,  Rodgers  (1968,  equation  (72))  has  shown 
that 


g 

AT  (W)  =  —  Ad  (w)  .  (9) 
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The  ALPS  displacment  amplitude  response  is  plotted  in  figure  44  equalized  to 
a  magnification  of  100,000  at  25  sec  period.  The  equivalent  tilt  response  is 
plotted  in  the  figure  to  the  same  scale  and  is  seen  to  reach  a  maximum  of  about 
30,000,000  for  a  rather  broad  peak  of  frequencies.  Rodgers  also  gives  an 
equation  for  the  output  of  a  seismograph  for  a  retrograde  Rayleigh  wave  as 


$  = 


(10) 


where 

$  =  actual  seismogram 
=  ideal  seismogram 
g  =  gravity  acceleration 
c  =  wave  phase  velocity 

6  =  particle  motion  ellipticity  parameter,  ratio  of  horizontal  semi¬ 
axis  to  vertical  semi-axis 

-f  =  period  of  wave. 

The  second  term  in  the  parenthesis  in  equation  10  is  the  tilt  term.  Values 
for  the  tilt  term  have  been  calculated  using  6  =  2/3,  which  would  apply  for 
an  infinite  half-space,  and  with  typical  values  for  Rayleigh  wave  phase 
velocity.  These  calculated  values  are  listed  in  table  3.  The  tilt  response 
g/c 

(1  +  - T)>  An  (dj)  ,  is  plotted  in  figure  44.  From  this  figure  it  can  be 

2tt6 

seen  that  the  addition  of  the  tilt  response  to  the  displacement  response  for 
a  Rayleigh  wave  is  very  small  compared  to  the  potential  response  that  a  horizontal 
seismometer  could  have  to  a  tilt  at  the  same  periods.  Ihe  difference  between 
the  Rayleigh  wave  response  and  the  tilt  response  will  be  referred  to  later. 

t 

il-6.2  Seismic  Spectral  Changes  with  Pressure  Activity  Changes 

Now  let's  consider  the  spectral  content  of  the  pressure  induced  noise,  then 
the  effect  on  the  horizontal  pendulum  of  pressure  induced  tilts  and  of  Rayleigh 
wave  induced  tilts  will  be  discussed.  Figure  45  shows  a  plot  of  the  S325L1 
spectra  under  low  and  high  wind  conditions  and  the  P325L  spectra  under  high 
wind  conditions.  As  the  wind  associated,  pressure  induced  noise  increases, 
it  increases  over  the  whole  spectra.  The  microseismic  peaks  visible  on  the 
S325L1  strain  during  the  low  wind  conditions  are  completely  swamped  with  the 
strains  from  the  high  wind  pressure  changes.  This  same  raising  and  lowering 
of  the  entire  spectrum  was  observed  by  Herrin  (1963,  a  personal  communication) 
with  a  vertical  seismograph  flat  to  velocity  from  1.25  sec  to  110  sec  used 
during  an  experiment  in  south  Texas.  At  low  wind  conditions  microseismic 
peaks  were  clearly  visible  at  2  Hz  and  6  sec.  As  the  wind  velocity  increased. 
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Figure  44. 


Horizontal  inertial  ALPS  seismograph  response 
and  Rayleigh  waves 


to  displacement,  tilt. 
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Table  3. 

Tilt  response  of  horizontal  inertial 
seismometers  to  Rayleigh  waves 

Period 

sec 

Phase  velocity 
km/sec 

&!s  T 

2 rB  T 

10 

3.0 

0.008 

15 

3.0 

0.012 

20 

3.36 

0.014 

24 

3.50 

0.016 

26 

3.56 

0.017 

30 

3.65 

0.019 

34 

3.  78 

0.021 

40 

3.90 

0.024 

45 

3.93 

0.027 

50 

3.97 

0.029 

55 

3.98 

0.032 

60 

3.99 

0.035 

70 

4.00 

0.041 

80 

4.03 

0.046 

90 

4.07 

0.052 

100 

4. 11 

0.057 

150 

4.28 

0.082 

200 

4.57 

0. 102 

2-0 

4.90 

f  t 

0.119 

300 

5.28 

0. 133 
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POWER  SPECTRAL  DENSITY  ((DISPLACEMENT  IN  M)2/H*  RELATIVE  TO  25  SEC.  SEE 


PERIOD  (sacondb) 


Figure  45.  Comparison  of  power  spectral  density  between  low  and 
high  wind  conditions 
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valleys  in  the  spectra  were  filled  in  until  finally  the  2  Hz  then  the  6  sec 
peaks  disappeared  into  the  pressure  induced  spectra.  Thus  the  pressure 
induced  noise  to  a  first  order  approximation  can  be  considered  to  be  a  white 
noise  source.  (The  ML2  spectra  in  figures  42  and  43  when  considered  along 
with  the  microbarograph  response  actually  indicate  a  flat  response  to  rate- 
of-change  of  pressure  between  about  50  sec  and  400  sec.)  The  output  of  any 
system  to  a  white  noise  source  is  the  system  response.  Therefore,  the  shape 
of  the  spectral  output  of  a  seismograph  can  give  strong  indications  as  to  the 
input  character. 

11.6.3  Non-Linear  Tilting 


Several  sets  of  data  demonstrating  the  tilt  response  have  already  been  pre¬ 
sented;  two  will  be  discussed  here.  The  spectra  of  the  inertial  P325L  during 
high  wind  conditions  are  compared  to  the  spectra  of  the  S325L1  strain  in 
figures  30  and  45.  At  periods  longer  than  64  sec,  the  S325L1  spectra  de¬ 
creases,  but  the  P325L  spectra  continues  to  increase.  This  divergence  in  the 
spectra  implies  a  non-linear  relationship  which  is  confirmed  in  the  coherence 
plotted  in  figure  30.  It  is  hypothesized  that  the  additional  output  on  the 
inertial  systen  is  a  symptom  of  rigid  body  tilting  with  possibly  relative 
movement  on  the  local  joint  and  fault  pattern.  Movement  on  local  joints  and 
faults  in  the  immediate  vicinity  of  instrument  systems  has  been  observed  by 
Harrison  (1969)  in  the  Poorman  Fault  near  Boulder,  Colorado.  Stewart  Smith  -• 
(1170,  personal  communication)  has  hypothesized  that  the  residual  strains 
observed  in  the  near  field  of  large  explosions  (and  possibly  in  the  far 
field  of  earthquakes)  are  actually  the  readjustment  of  the  strain  field  in 
the  vicinity  (10's  m  to  10's  km)  of  the  seismometers. 

For  teleseismic  earthquakes,  this  readjustment  is  triggered  by  the  passage  of 
the  Rayleigh  wave  which  is  generally  the  largest  amplitude  elastic  wave. 

Both  Harrison's  and  Smith's  observations  apply  to  a  step  readjustment.  The  . 
hypothesis  presented  here  is  that  the  readjustments  along  local  joints  and 
faults  is  a  dynamic  phenomena  occurring  all  the  time, in  much  smaller  ampli¬ 
tudes  (on  the  order  of  0.1  Angstrom)  than  have  previously  been  resolved.  It 
is  felt  that  this  phenomena  is  not  unique  to  QC-AZ  but  that  it  can  be 
observed  at  any  high  sensitivity  seismic  observatory. 

The  LR1  Rayleigh  wave  spectra  in  figlire  40* also  substantiates  the  excess  tilt 
response  of  the  horizontal  inertial  P325L  oveir  the  response  of  the  strain 
S325L1.  If  the  S325L1  spectra  is  taken  as  representative r,of  the  true  Rayleigh 
wave  motion,  then  table  3  and  figure  44  indicate  that  at  100  sec  there  should 
only  be  about  5.7  percent  greater  amplitude ~on  the  pendulum  from  the  Rayleigh 
wave  tilt.  The  assumption  that  the  S325L1  spectra  represents  the  true  Ray¬ 
leigh  wave  motion  is  supported  by  the  cross  spectra  in  figure  40  which  is 
almost  identical  to  ;the  S3251M  spectra.  The  observed  output  in  power  spectral 
density  at  100”  sec  of  the  P325L  is  a  factor  of  about  30  larger  than  the  power 
of  the  S325L1.  The  ratio  of  the  observed  30:1  difference  is  so  much  greater 
than  the  1.057:1  difference  predicted  by  the  theory,  that  it  has  to  be  related 
to  a  phenomenon  other  than  an  elastic  response  to  the  travelling  wave.  Even  < 

if  the  S325L1  spectra  is  nqrmalized  upward  by  the  3.5  power  ratio,  the  tilt 
output  is  still  in  excess  over  the  Rayleigh  wave  tilt. 


In  sjmmary,  the  authors  contention  is  that  this  other  phenomenon  is  a  rigid 
block  motion  on  the  order  of  0.1  A  in  amplitude  on  the  local  joint  pattern 
and  that  this  motion  is  apparently  accompanied  by  a  frictional  damping  action 
It  is  expected  that  this  same  action  will  be  observed  at  any  high-sensitivity' 
seismic  observatory.  6  »  3 


12.  SPECIAL  REPORT 
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\  ,  , 

A  special  technical  report  (Fix  and  Sherwin,  1970)  has  been  submitted  to 
summarize  the  instrumentation  complex  objectives,  design,  installation,  and 
preliminary  results.  Authority  for  publication  has  been  received  according 
H  ^1FOf  security  review  proceduie.  This  paper  has  been  submitted 

to  the  Bulletin  of  the  Seisirological  Society  of  America  for  consideration  of 
publication. 
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INSTRUMENT  NOMENCLATURE,  RECORDER  CHANNEL  ALLOCATION.  AND  POLARITY 


This  appendix  lists  the  instrument  nomenclature,  the  normal  recorder  channel 
assignments,  and  the  polarity  convention  for  the  Queen  Creek  Seismo logical 
Station  (QC-AZ) .  6 


1.  NOMENCLATURE 

The  nomenclature  used  for  the  instrument  system  identifiers  is 

Category  of  data  Identifier 

Time:  BCD  =  binary  coded  decimal  station  time 

WWV  =  WWV  time 


Compensation:  Comp  =  playback  compensation  channel 


Data: 


first  symbol:  S  =  strain 

P  =  pendulum 
M  *  microbarograph* 
W  =  wind 


second  symbol*:  Z  =  vertical 

325  =  325  degree  azimuth 
55  =  55  degree  azimuth 
V  =  (wind)  velocity 
D  -=  (wind)  direction 


third  symbol:  L  =  Advanced  Long-Period  System  (ALPS) 

response 

U  =  Ultra-Long-Period  (ULP)  response 
X  =  Extended-Long-^eriod  response 
B  =  Broad-Band  (BB)  response 
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c- 


S  =  Short-Period  (SP)  response 

fourth  symbol:  L  =  Low  gain 

1,2,3  =  numbers  assigned  to  similar 
instruments 


*Microbarograph  channel  identifiers  do  not  have  a  symbol 
in  the  second  symbol  category. 

flag  =  identifier  for  clipping  at  Ithaco 
amplifier  output 


2.  RECORDER  CHANNELS 


2 . 1  TAPE  RECORDERS 


Channel  No. 

No.  1 

Recorder 

No.  2 

No.  3 

1 

BCD 

BCD 

BCD 

2 

SZLL 

WV2 

MU 

3 

P325LL 

SZB 

SZUL 

4 

P55L 

P55S 

P55X 

5 

S325LL 

S55B 

S325UL 

6 

S55L 

S55S 

S55U 

7 

Comp 

Comp 

Comp 

8 

S325L 

S325S 

S325U 

9 

SZL 

SZS 

SZU 

10 

P325L 

P325S 

P325X 

11 

PZL 

PZS 

PZX 

12 

S55LL 

WD2 

S55UL 

13 

PZLL 

S325B 

ML1 

14 

P55LL 

WWV 

ML2 

TR  70-29,  app. 

II*  IK  i  fill  ***fcA* 


-2- 


a 

1 


2.2  16  mm  FILM  RECORDERS 


OJ 


Channel  No 

1 

2 

* 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 


No.  1 


Re corder 


No.  2" 


/ 


FLAG 

WWV 

SZL 

SZLL 

S325L 

S325LL 

S55L 

S55LL 

PZL 

P32SL 

P55L 

PZX 

P325X 

P55X 

ML1 

PZS 

P325S 

P5SS 

WV 

WD 


Wh'V 

S325L-S55L 

SZL+S525L+S55L 

S325L+S55 L 

SZL 

PZL 

S325L+P325L 

S325L 

P325L 

S32SL-P325L 

S55L+P55L 

S55L 

P55L 

S55L-P55L 

ML1 

off 

only 

16 

channels 


2.3  35  mm  FILM  RECORDER 


Channel  No. 

1 

2 

3 

4 


szu 

S55U 

S325U 

MU 


3.  POLARITY 


The  polarity  convention  used  on  QC-AZ  recorders  is  as  follows: 


Strain: 


Compression  is  up  on  film  recorders  and 
positive  voltage  into  the  magnetic-tape 

T*P  rno  rc  * 
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Inertial: 


Displacement  toward  up,  toward  525 
deg  azimuth,  and  toward  55  deg  azimuth 
are  up  on  film  recorders  and  positive 
voltage  into  the  magnetic-tape  recorders. 


Electromagnetic  calibration 

coils:  Connected  so  that  a  dc  current  turned 

,  on  produces  a  positive  output. 


-r 


i 


C 


■v 
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ABSTRACT 


three  strain 

m*n*'  A  -eond  Develacorder  LZ  “«*  «• 

seismoyaphs.  The  second  325  deg  azimuth  strain  seismometer  w’iuJ"  J0"3  k®3"1*  °n'line  by  comb,n,n9  «Rin  and  inertial 
- coils  CO  8,.  LP  *anU  ZLZZZZZTT*  “  Cmttua  ‘  ,OTiC*1  MNn«. 

filter  on  the  input  to  the  strain  preamplifier  has  limited  amni  f-  °R  l9Uratlon  to  match  ^  strain  response,  A  low-pass 
sea.  has  been  improved  to  a  4  ^  **  ++  ""  °f  ™e  rfine 

»wre  investigated.  The  most  practical  is  to  record  the  outout  fr  °  s*31™-"9  svstem  transfer  ^ 00 1  ions  from  special  test  data 
determine  the  transfer  function  by  numerical  smoothina  and  d  ff"  * '"’pulse  of  current  into  ^  calibration  coil  and 

provide  enhancement  of  signals,  wave  type  identification  and  'hj _ ia*inB-  agnetic  tape  playback  and  on-line  recordings 

teleseismic  earthquakes.  Enhanced  short-period  recording  *ima  fr<>m  mterferrin9  events  and  from  low-level 

Colorado,  to  QC-AZ  are  sign.^ft.y  59  ^  Ra" 

mi,  =  4.1  at  62  deg.  mb  =  4.0 at  25 deg^Jd  mh  =^2^M00deo  Ex*"?les  of  smaM  waves detected  are: 

parisons  between  parallel  strain  seismometers  durini  low  and  hir*  densitY'  coherence,  and  phase  angle  corn- 

instruments  are  operating  proper! /and  that  calibrations  are  Z  C°^d'tlons  and  earthquake  signals  demonstrate  that  the 

is  reviewed.  Theoretical  °f  **  eVa'U3ti°n  °f ' -nd  related  LP  noise 

50  sec.  a  horizontal  inertial  seismometer  operating  at  100K  («  25  ^iUSwffdT'  Y'T  T*  fluctuation-  At 

Rayleigh  wives  of  46,300  and  for  tilts  of  27  900  000  Ahv^tw  "^cations  for  d.splacements  of  45,000,  for 

non-linear  rigid  block  tiltjng  occurs  as  a  result  of  Dressure  l!^d  Ha  ^  31 3!'  hard  r<Xk  se'smo9raPh  stations 

is:  (1)  There  is  a  slight  relative  amplitude  shift  with  frequencv  and  a  90^  wav“- Evidence  *at  "PPorts  this  hypothesis 
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